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SUMMARY 


The use of eigenspace assignment techniques to synthesize flight control 
systems for flexible aircraft is explored. These eigenspace assignment tech- 
niques are used to achieve a specified desired eigenspace, chosen to yield desir- 
able system impulse residue magnitudes for selected system responses. Two 
eigenspace assignment techniques are investigated. The first is a technique for 
directly determining constant measurement feedback gains that will yield a 
closed-loop system eigenspace “close” to a desired eigenspace. The second 
technique is a method for selecting quadratic weighting matrices in a linear 
quadratic control synthesis that will asymptotically yield the closed-loop 
achievable eigenspace. Finally, the possibility of using either of these tech- 
niques with state estimation is explored. 

Application of the methods to synthesize integrated flight- control and 
structural-mode-control laws for a large flexible aircraft is demonstrated and 
results discussed. Eigenspace selection criteria based upon the design goals are 
discussed, and for the study case it would appear that a desirable eigenspace 
can be obtained. In addition, the importance of state-space selection is noted 
along with problems with reduced-order measurement feedback. Since the 
full-state control laws may be implemented with dynamic compensation (state 
estimation), the use of reduced-order measurement feedback is less desirable. 
This is especially true since no change in the transient response from the pilot’s 
input results if state estimation is used appropriately. The potential is also 
noted for high actuator bandwidth requirements if the linear quadratic syn- 
thesis approach is utilized. Even with the actuator pole location selected, a 
problem with unmodeled modes is noted due to high bandwidth. 

Some suggestions for future research include investigating how to choose 
an eigenspace that will achieve certain desired dynamics and stability robust- 
ness, determining how the choice of measurements effects synthesis results, and 
exploring how the phase relationships between desired eigenvector elements 
effects the synthesis results. 
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CHAPTER I 
INTRODUCTION 


Many aircraft control systems in use today are synthesized by treating the 
aircraft as if it were a rigid-body. This is appropriate when the aircraft are 
rigid enough so that the structural flexibility has little effect on the rigid-body 
dynamics. Because of the rapidly growing use of composites, the need for 
greater payload capability and increased range, future aircraft will become 
lighter and more flexible. Therefore, the aircraft’s aeroelastic modes will begin 
to play a greater role in its dynamics. 

These aeroelastic modes pose special problems with both the analysis of 
the flexible dynamics and in synthesizing control systems to achieve desired 
dynamics. Work done by D.K. Schmidt and M.R. Waszak [1,2,3] addressed the 
first problem, analyzing the dynamics of aircraft in which structural flexibility 
is significant. This report will build upon the work of Schmidt and Waszak 
and address the problem of how to synthesize control systems to achieve 
desired dynamics in flexible aircraft. In particular, these control systems will 
be synthesized using eigenspace assignment techniques [4, 5, 6, 7, 8]. This work 
will first explore two eigenspace assignment techniques then examine how these 
techniques can be used with state estimation. 

This report is presented in six chapters. A brief summary of Schmidt and 
Waszak’s work and how this work naturally leads to an eigenspace assignment 
control synthesis technique is presented in Chapter II. Chapter HI presents 
Direct Eigenspace Assignment (DEA). DEA is a control synthesis technique for 
directly determining measurement feedback control gains that will yield a 
closed-loop system eigenspace as close as possible to a desired eigenspace 
[5, 6, 7, 8]. In Chapter IV, Linear Quadratic Eigenspace Assignment (LQEA) is 
presented. This method is based upon the work of C.A. Harvey and G. Stein 
[9,10]. LQEA is a method for selecting quadratic weighting matrices that will 
asymptotically yield a closed- loop achievable eigenspace in a linear quadratic 
control synthesis. Chapter V explores the advantages of using either of these 
eigenspace assignment techniques with state estimation. In conclusion, Chapter 
VI presents a summary of the results and conclusions based on the results. 
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CHAPTER H 
BACKGROUND 


A control synthesis technique can be viewed as a mapping from a given set 
of poor dynamics into a set of desirable dynamics (See Figure 2.1). When con- 
sidering a rigid aircraft, these desirable dynamics are well defined. They are 
stated in terms of handling qualities specifications, time response requirements, 
etc. These desirable dynamics can usually be expressed in a way that is com- 
patible with a popular control synthesis technique, such as Root Locus, Fre- 
quency Analysis methods (Bode, Nyquist, etc), or Linear Quadratic optimiza- 
tion methods. Therefore, using one of these techniques a designer can usually 
synthesize a control system for a rigid aircraft that will achieve some desired 
dynamics in the closed-loop system. 

In reality, aircraft are not rigid. The term “relatively rigid” only refers to 
aircraft in which the vibrational frequencies of the aircraft’s aeroelastic modes 
are large compared to the frequencies of the rigid-body modes. As the frequen- 
cies of the aeroelastic modes becomes lower, these aeroelastic modes begin to 
significantly affect the aircraft’s dynamics. These structural modes pose special 
problems, not only with describing the flexible dynamics, but also in synthesiz- 
ing control designs to achieve desired dynamics. 

One control synthesis approach that is too often used is just ignoring the 
flexible dynamics. But work done by M.G. Gilbert, D.K. Schmidt, and T.A. 
Weisshaar [11] has shown that ignoring the structural dynamics of these air- 
craft in flight control designs can lead to very poor if not disastrous results. 


Modal Analysis 

In work on the analysis of flexible aircraft dynamics, D.K. Schmidt and 
M.R. Waszak [1,2,3] present one solution to the problem of describing the 
effects of flexible dynamics. Schmidt and Waszak present a modal analysis 
technique that is especially useful for analyzing the dynamics of flexible vehi- 
cles. This technique involves calculating both the eigenvalues and impulse resi- 
dues of the system. By comparing the magnitudes of the system’s impulse 





Desired 

Configuration 


Figure 2.1 

Control Synthesis Technique Mapping 
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residues for various outputs, this analysis shows how much the elastic modes 
contribute to the system’s dynamics. The residues of a system are related to 
its dynamics in the following way. 

Given the system 

x = Ax + Bu (system dynamics) 


y = Cx (system responses) 

where x e R n , u e R m , and y . e R p . Assume all the eigenvalues of A are dis- 
tinct. Transform this system to modal coordinates using the transformation 

x = Tq 


where T is the system modal matrix (a matrix of the eigenvectors of A). One 
obtains 

4 = T^ATq + T -1 Bn 


= A 4 + Bu 


where A is a diagonal matrix of the eigenvalues of A, and B is the system con- 
trollability matrix. Taking the Laplace transform of this equation (assuming 
q(0) = 0) yields 

sq(s) = Aq(s) + Bu(s) 


a(s) = [sl„ - A]-' Bu(s) 


The system input-output transfer function relationship is given by 

S (s) = CTq(s) = C[sl„ - A]" 1 Bu(s) 


Looking at the ith output, jth input 


Yi(s) 

Uj(s) 


= Cj [sl n - A]" 1 b: 
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_ y ^i,k^k,j _ y, C i t/ k 1 k b j 

k=l S " X k k = l S " X k 

where 

Cj = i th row of C 
Cj = i th row of C 
bj = j th column of B 
bj = j th column of B 
c ik = (i,k) element of C 
b k j = (k,j) element of B 
X k = k fch eigenvalue of A 

Uy = k th column of T (k th right eigenvector of A) 
l k = k th row of T _1 (k th left eigenvector of A) 

For an impulse in the j th input, the i th output is given by 

Uj(t) = impulse — ► Uj(s) = 1 


yi(s) 


y 6 i,^kj _ “ 

k=! S_X k k=l S ~ X k 


n i k 
k=i s “ X k 

The system output time response due to an impulsive input can be found 
by taking the inverse Laplace transform of this equation. This yields 

Yi (t) = t Ri,j,k e x “ 

k=l 

In this equation R;j jk is the impulse residue for output i, associated with 
eigenvalue k, and due to an impulse in input j. 
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Example 2.1 

This example demonstrates the role impulse residues play in a flexible 
aircraft’s pitch rate dynamic response due to an impulsive elevator input. For 
this example, let y;(t) = 0(t) and Uj(t) = <5 e (t) 

m = e R i eX, ‘ 

i=l 


= R P h e Xpht + R ph + R sp e^ + R sp e^ 

+ E ( R 6 e X( ‘ + ) 


= 2 1 R Ph | e <Tpht cos ("pht + $ P h) + 


where 


2 1 R sp | e ff ‘ pt cos (w sp t + $ sp ) + 

E 2 1 R fc| e cos K* + $&) 

i=l 


Xj = -o - ; + jw; , 4>j = arctan ( Im(R;)/Re(Rj) ), 
| Rj| = (Re 2 (R;) + Im 2 (Rj))R 2 


k = number of elastic modes 
sp = short period mode 
ph = phugoid mode 
= i th elastic mode 
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As this example shows, a system’s dynamics are dependent on both its 
eigenvalues and its residues. The eigenvalues determine the natural frequency 
and damping of each mode. The residues are an indicator of how much each 
mode of the system contributes to a given output. 

Schmidt and Waszak’s research involved analyzing a family of 
configurations of large flexible aircraft similar in geometry to the Rockwell B-l 
(see Figure 2.2). The structural rigidity of these configurations was varied by 
parametrically varying the invacuo-structural vibration frequencies of the two 
lowest frequency aeroelastic modes of the vehicle. A summary of the analysis 
of two of the configurations will be presented to illustrate the method. 


Example 2.2 

Configuration One (Baseline) 

This configuration is representative of an unmodified B-l. This 
configuration received good subjective ratings of the handling qualities t in a 
tracking task with pilot comments of “Good; no problem”. The eigenvalues of 
this configuration are: 

X sp = - 1.49 ± j 2.37 = (w n = 2.799 (rad/s) , f= 0.532) 


X ph = - 0.0015 ± j 0.0672 = (« n = 0.067 (rad/s) , j= 0.022) 


X €l = - 0.657 ± j 13.29 = (w n = 13.31 (rad/s) , f= 0.049) 


X| 2 = - 0.46 ± j 21.35 = (w n = 21.35 (rad/s) , 0.022) 


The impulse residue magnitudes (expressed as a percentage of the total 
response) of the rigid-body pitch and pitch rate outputs can be found in Figure 


2.3. 


f The Cooper Harper rating scale is a subjective rating (1 being best and 10 being worst) 
used to describe vehicle handling qualities in various tasks [12]. 
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Figure 2.3 

Impulse Residue Magnitudes - Configuration One 
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Figure 2.4 

Impulse Residue Magnitudes - Configuration Two 
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Configuration Two 

This configuration differs from Configuration One in that it has a lower 
first aeroelastic mode invacuo frequency. Therefore, this configuration is more 
flexible than Configuration One. The pilot ratings and tracking scores this 
configuration recieved (in the tracking task) were worse than those of 
Configuration One. The pilot comments were “little oscillation; slight control 
response lag”. The eigenvalues of this configuration are: 

A sp = - 1.348 ± j 2.193 = (u; n = 2.574 (rad/s) , 0.524) 


\ ph = - 5.71 x 10~ 5 ± j 0.053 = (w n = 0.053 (rad/s) , £= 0.001) 


A fi = - 0.726 ± j 8.757 = {oj n = 8.787 (rad/s) , £= 0.083) 


X^ 2 = - 0.456 ± j 21.351 = (ou n = 21.35 (rad/s) , ?= 0.021) 


The impulse residue magnitudes of the rigid-body pitch and pitch rate outputs 
can be found in Figure 2.4. By comparing the impulse residue magnitudes of 
these two configurations one can see that Configuration Two, which received a 
lower pilot rating, has a greater contribution from the first elastic mode in 
these outputs. For a more detailed explanation of this method the reader is 
referred to Reference [3]. 


Schmidt and Waszak summarize their modal analysis research results by 
stating: “The results of modal analysis indicate when the magnitudes of the 
modal impulse residues of the aeroelastic modes become large compared to the 
residue magnitudes of the rigid-body modes for important outputs, the dynam- 
ics can change significantly and in such a way that the handling qualities of the 
vehicle may be degraded. 


f Reference [3], p.85. 
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Control Synthesis 

This analysis leads to one possible approach to the flexible aircraft control 
synthesis problem. This approach would be to synthesize a control system that 
would control both the eigenvalue locations and impulse residue magnitudes 
associated with undesirable modes in certain outputs. For example, reducing 
the magnitude of impulse residues associated with aeroelastic modes in rigid- 
body outputs would make the flexible aircraft behave more like a rigid aircraft. 
This type of control strategy would reduce the undesirable effects of structural 
flexibility in the aircraft’s dynamics. 

This motivates interest in a control synthesis technique that will modify 
both the system’s eigenvalues and its impulse residue magnitudes. But, as was 
shown, impulse residues are a function of the system’s eigenvectors. This 
naturally leads to a control synthesis technique that involves achieving some 
desired eigenspace in the closed-loop system. This type of control synthesis 
technique will be referred to as an eigenspace assignment technique. 

In the following chapters, two eigenspace assignment techniques will be 
presented. This report will investigate each of these with regard to level of 
design complexity, the utility afforded to the designer, and how well each of 
these techniques enhances the vehicle’s dynamics. Each technique will be 
evaluated, demonstrated in examples, and compared to see how well it meets 
each of these control synthesis goals. 
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CHAPTER HI 

DIRECT EIGENSPACE ASSIGNMENT 


In this chapter, a control synthesis technique for directly determining 
measurement feedback control gains that will yield an achievable eigenspace in 
the closed-loop will be developed. This technique will be referred to as the 
Direct Eigenspace Assignment (DEA) technique. 

In papers on eigenspace assignment S. Srinathkumar [5,6]; A.N. Andry, 
E.Y. Shapiro, and J.C. Chung [7]; and T.B. Cunningham [8]; after B.C. Moore 
[4], have explored the limits of the achievable eigenspace for a system with 
measurement feedback. For a system that is observable and controllable and 
has n states, m controls, and / measurements; one can exactly place / eigen- 
values and m elements of their associated eigenvectors in the closed-loop sys- 
tem t. If it is desired to place q elements (m < q < n) of the / eigenvectors 
associated with the / desired eigenvalues, the “best” that can be achieved is a 
least squares fit to the desired eigenvectors. 

The following section will present the development of the DEA synthesis 
technique. Later sections will present a discussion of the achievable eigenspace 
and three examples that illustrate some of the strong and weak points of DEA. 


DEA Technique Development 

Given a system 

x = Ax + Bn (system dynamics) (3.1a) 

y = Cx (system responses) (3.1b) 


t This assumes l > m. For a general statement and proof of this property the reader is 
referred to Reference J5]. 
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2 . = Mx. 

(system measurements) 

(3.1c) 

Mc= G 1 

(feedback control law) 

(3.1d) 

U = Hc+ Up 

(total control input) 

(3.1e) 


u p = pilot' s inputs 

where x e R n , u c R m , y e /? p , and z e R l . One can express the spectral 
decomposition of the open-loop system, A, by 

A .&= X^. i=l,...,n (3.2) 


where 

X 0 . = i th open-loop eigenvalue 

— i th open-loop eigenvector 

The system augmented with the control law is 

i = (A + BGM)x + Bu p (3.3) 

and therefore the spectral decomposition of the augmented (closed-loop) system 
is given by 

(A + BGM)^ = X cA i=l,...,n (3.4) 

Arranging this, one can obtain 

BGMi/ c = X^- A^. (3.5) 

In Equation (3.5) let 

wj = GM^ (3.6) 


and then solve for v... One obtains 

— -c, 
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Bw 5 = X Cj ^- A^. 

= A)^. 

and finally 

2c f = (^cJn- A) _1 B5Ki (3.7) 


= eigenvector achievable in the closed-loop 
system for a given X c . and Wj 


This equation describes the achievable i th eigenvector of the closed- k 
system as a function of the closed-loop eigenvalue, X c ., and Wj. 

If one could calculate a value for w; that would make the achievable eigen- 
vector, v a ., as close as possible to some desired eigenvector, it could be used 
to determine a gain matrix that would yield this eigenspace. This is done by 
defining a cost function associated with the i th mode of the system (assume for 
now that X c ., f/ d ., and are real) 

J i= - tear kd^QdKr i=1 > •>* te-s) 


where 

l/ a . = i th achievable eigenvector associated 
with eigenvalue \ 

v_ A . = i th desired eigenvector 

Qd; — i th n-by-n symmetric positive semi-definite 
weighting matrix on eigenvector elements 

This cost function represents the error between the achievable eigenvector and 
the desired eigenvector weighted by the matrix Q d . 
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By substituting the equation for the achievable closed-loop system eigen 
vector, Equation (3.7), into the cost function, Equation (3.8), one obtains 

h = J tea, - ^Qditea, “ I'd,) 

= j- ((( x i I n - A] _1 Bwj - ^.) T Q d .([X ; I n - A] _1 Bwi - ^.)) 

= j (wi T B T (XjI n - A]- T Q di [XiI n - A] -1 Bffij 
“ ^d^QdJXiln - A] _ 1 Bw| - w i T B T [X i I n - A]~ T Q dj i/ d . 

+ lid, T QdA ) 


To find the W; that will minimize J;, take the gradient of J; with respect to 
Wj and set it equal to zero. 



= - Ar T Q dl |X,I„ - A]-‘B 


£d, T Qd,M» ~ A]-'B 


(3.0) 


Now solve Equation (3.9) for W; T 



- i^d i T QdI j ili , i T Qd, Li] 1 


(3.10) 


where 
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Lj = (X d .I n - A) -1 B 
X d . = i th desired eigenvalue 
£ d . = i th desired eigenvector 
Q dj = i th n-by-n symmetric positive semi-definite 
weighting matrix on eigenvector elements 

Note that in this development X d . cannot belong to the spectrum of A. 

By examing Equation (3.6), one can see that given values of w; and i/ a . for 
i = l, . . . , / one can determine the feedback gain matrix, G. This is done by 
using Equation (3.10) to calculate Wj’s for i=l,...,/. Concatenate the individual 
wfs columnwise to form W. 

W = [ wj | w 2 | • • • | w, ) 

Use Equation (3.7) to calculate the closed-loop system achievable eigenvalues 
j/ a .’s for i=l, ...,/. Concatenate the individual, i/ a .’s, columnwise to form V. 

V = U ai | ^ | • • • | £a, } 


One can now solve for the gain matrix that will yield the l desired eigen- 
values and / achievable eigenvectors in the closed-loop system by applying 
Equation (3.6). 

[ Hi | • • • | YLi ] = GM[ u %i | • • • | v Zi ] 


W = GMV 


Solving for G, one obtains 

G = W(MV] _1 


(3.11) 


The augmented system is given by 
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i = A aug x. + Ba p (3.12) 

where 

A aug = A + BGM (3.13) 


The augmented system transfer functions are obtained from 

sk(s) = A aug x.(s) + BHp(s) 


(sl n ~ A aug )x(s) = Bu p (s) 


x(s) = (sl n - A aug ) *6^(8) 


The system responses are given by 


y(s) = Cx(s) 


(3.14) 


therefore, 

y(s) = C (sl n - A aug ) _1 B Up(s) (3.15) 


The solution method is similar for the case of control feed through in the 
measurements (i.e. z = Mx + Nu) and for the case of complex eigenvalues and 
eigenvectors. The DEA control synthesis technique developments for both of 
these cases are given in Appendix A.l. 


Achievable Eigenvectors 

As was shown in the DEA development, Equation (3.7) 

14*= Af'Bw , (3.7) 

describes the achievable eigenspace of the system. In this equation, 
(\ c .I n - A) -1 B is an n-by-m matrix and w ; is an m-by-1 vector. By examing 
this equation, one can see that the number of control variables (m) determines 
the dimension of the subspace in which the achievable eigenvectors must reside. 



18 


If the desired eigenvector, u j., cannot be achieved, the closest one can come to 
it (the achievable eigenvector) is a projection of onto the subspace spanned 
by the columns of (Xjl n — A) -1 B (see Figure 3.1). 

If one desires to exactly place a larger number of eigenvector elements, the 
rank of B must be increased (i.e. one must have more controls). If one is 
interested in having control over a larger number of modes (exactly placing 
eigenvalues and a least squares fit for the eigenvectors), the rank of M must be 
increased (i.e. more measurements). 


System Scaling 

In general, the state vector must be transformed and the units changed in 
the following way to aid in the physical interpretation of the eigenvector ele- 
ments. Let x s be the new state definition, where 

x s = Tx 


then 


x g = TAr’xs + T Bu 


where 


i=cr'L 


i\i r'i. 


T — diagonal matrix of scaling elements 


For the large flexible aircraft examples presented in the following section, T is 
chosen to be 

T = diag(1.0, 1.0, 1/U 0 , 1.0, </>' , , <^' 2 , <j>> x , <j>< 2 ) 


with 


U 0 = forward flight speed (ft/sec) 
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h 



Figure 3.1 

Achievable Subspace 
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<j >' ; = mode slope of the i th elastic mode 
at the cockpit 

The scaling matrix used in the examples is given in Appendix A.2. All of the 
example results are presented in the scaled units. The reader should note that 
a system’s responses, eigenvalues, and impulse residues are independent of sys- 
tem transformations, but the eigenvectors are not. Therefore, the choice of the 
system scaling will effect the choice of the desired eigenvectors and the syn- 
thesis results. 


DEA Examples 

Three examples will be presented to demonstrate the method. All three 
will use the same model, the longitudinal dynamics of a large flexible aircraft 
(Schmidt and Waszak’s Configuration Two). The model is based on a steady- 
state flight condition of forward cruise speed equaling 949 feet/second at sea 
level. It includes the four standard rigid-body degrees of freedom and the two 
lowest frequency symmetric aeroelastic modes of the vehicle. These two aeroe- 
lastic modes are the first fuselage bending mode and the second fuselage bend- 
ing mode. 


System 

The unsealed model is as follows: 


X = Ax + Bn (system dynamics) 


y. = Cx + Du. (system responses) 


z = Mx + Nu (system measurements) 


li = He + Up 


(total control input) 
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u c = G 2 . (feedback control law) 

Up = pilot' s input 

States 

x T = ( a , 0 r , u f , 0 r , ^ , £1 , £2 ) 

where 

a = angle of attack (radians) 

9 r — rigid-body pitch rate (radians/sec) 
u f = forward velocity (ft/sec) 

6 X — rigid-body pitch angle (radians) 

= mode one generalized deflection (dimensionless) 
£2 = mode two generalized deflection (dimensionless) 
= mode one generalized deflection rate (1/sec) 

£2 = mode two generalized deflection rate (1/sec) 

Controls 

Hc T = ( . 5 cv ) 

llp T = ( <5 e , 0.0 ) 

where 


<5 e = elevator deflection (radians) 
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5 CV = forward control vane deflection (radians) 


Outputs 

X T = ( , 0 r , 9 r , 7 , 9 t , 9 t , n Zp ) 

where 

Uf = forward velocity (ft /sec) 

0 T = rigid-body pitch angle (radians) 

0 t = rigid-body pitch rate (radians/sec) 

7 = flight path angle (radians) 

0 t = total pitch angle (radians) 

0 t = total pitch rate (radians/sec) 

n Zp = plunge acceleration at the cockpit (g 1 s) 

The total pitch angle is the pitch attitude measured at the cockpit. It is 
the rigid-body pitch attitude plus contributions of the local structural 
deflections (Figure 3.2). The total pitch angle is defined by the equation 

», = «,- E«'i(U 

i=l 

where 

k = number of aeroelastic modes 

l x = distance between the c.g. and station x (ft) 

(j )' ; = mode slope of the i th elastic mode (ft/ft) 

= generalized coordinate of the i th elastic mode 


The plunge acceleration at location x, n Zx , is defined by 
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»* = ■—*■ I u.t + V» t " EfA(><) I 

S i=l 

where 

g = gravitational acceleration (ft/sec 2 ) 

U 0 = cruise velocity (ft /sec) 

l x = distance between the c.g. and station x (ft) 

<f > ; = mode shape of the i th elastic mode (ft) 

The flight path angle, 7 , is given by 

1 = 0 r - a 


Measurements 

The measurements considered for feedback are 

1 T = ( . K > n 2p , n Zt ) 

where 

9 t = total pitch angle measured at the cockpit (radians) 

0 t — total pitch rate measured at the cockpit (radians/sec) 

n Zp = plunge acceleration at the cockpit (g / s) 

n Zt = plunge acceleration at an aft fuselage station (g' s) 


For the examples, the aft fuselage station was chosen to be located at 
l x = -28.23 ft. The open-loop system eigenvalues, the eigenvalues of A, are: 

\ sp - - 1.348 ± j 2.193 = (w n = 2.574 (rad/s) , 0.524) 
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X ph = - 5.71 x 10 5 ± j 0.053 = (w n = 0.053 (rad/s) , 0.001) 

X^= - 0.726 ± j 8.757 = (w n = 8.787 (rad/s) , 0.083) 

X <2 = - 0.456 ± j 21.351 = (w n = 21.35 (rad/s) , $= 0.021) 


The open-loop system eigenvectors are given in Table 3.1. In this table, the 
eigenvectors have been scaled such that the first element is unity. The open- 
loop system A, B, C, D, M, and N matrices are given in Appendix A.2. The 
open-loop system impulse residue magnitudes for each output due to each mode 
are shown in Figure 3.3. The residue magnitudes for each response are normal- 
ized for plotting such that the sum of the residue magnitudes for each response 
is unity. By examining the relative impulse residue magnitudes one can note 
several things. The first elastic mode contributes significantly to the rigid-body 
pitch rate, 0 t . The first elastic mode also contributes to the rigid-body pitch 
angle, 0 r . The elastic modes contribute significantly to the total pitch rate, 0 t , 
total pitch angle, 0 t , and plunge acceleration at the cockpit, n Zp , outputs. 

The open-loop system time responses due to a unit step in u p , with initial 
conditions x(0) = Q, are shown in Figure 3.4. The rigid-body pitch rate time 
response also shows the first elastic modes contribution to this output. The 
high frequency ripple is due to the first elastic mode. 


Example 3.1 : Four Measurements 

This control synthesis has two design objectives - to improve the aircraft’s 
rigid-body dynamics and to reduce the undesirable effects of flexibility on the 
rigid-body dynamics. For this example dim(x.)— n=8, dim(n)=m=2, and 
dim(z)=/ : =4; therefore using DEA one can modify four eigenvalue/eigenvector 
pairs. 


Specification of Synthesis Parameters 

Desired Eigenvalues. For this control synthesis one may specify four real 
desired eigenvalues or two complex conjugate eigenvalue pairs. By examining 
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Table 3.1 

Open-loop Eigenvectors 


CONFIGURATION TUO : OPEN LOOP 


EIGENUALUES 
1 : - 1.3479 

2 J - 1.3479 

3 : -.0001 

4 : -.0001 


+ J 2.1929 

- J 2.1929 

+ J .0527 

- J . 0527 


EIGENUECTORS: (MAGNITUDE. PHASE ( DEGREES ) ) 


1 


( 

l.OOOOE+OO* 

0) 

( 

1.0000E+00, 

( 

2.2949E+00, 

-2.6416E+02) 

( 

2.2949E+00, 

c 

2. 1545fc“02> 

-3. 1488E+02) 

C 

2.1545E-02, 

c 

8.9153E-01, 

-2.5741E+01) 

( 

8.9153E-01 * 

( 

2 . 8846E-0 1 ♦ 

“3.3438E+02) 

C 

2.8846E-G1 9 

( 

5.1928E-02, 

-1.7132E+02) 

( 

5. 1928E-02t 

c 

7. 4252E-01 » 

-2. 1280E+02) 

( 

7.4252E-01, 

( 

1.33BBE-01, 

-4.9746E+01) 

( 

1.3366E-01, 

c 

1 . 0000E+00f 

3 

0) 

( 

1.0000E+00, 

( 

1 .797BE-01* 

-1.9265E+02) 

( 

1.7976E-01, 

( 

2.1292E+00, 

-1.7933E+02) 

( 

2.1292E+00, 

( 

3.4097E+00, 

-2.8271E+02) 

( 

3.4097E+00* 

( 

2.527SE-01* 

3.6763E-01) 

( 

2*5276E-01» 

c 

5.2150E-02, 

~1 • 7985E+02) 

c 

5.2150E-02* 

( 

1.3326E-02, 

“2.6957E+02) 

c 

1.3326E-02, 

c 

2.7494E-03, 

-8.9783E+01 ) 

( 

2.7494E-03, 


2 

0) 

2.G416E+02) 
3. 1488E+02) 
2. 5741E+01 ) 
3.3438E+02) 
1.7132E+02) 
2. 1280E+02) 
4.974GE+01 ) 


4 

0 ) 

1.92G5E+02) 
1 . 7S33E+02) 
2. 8271E+02) 
3.G7G3E-01) 
1.7985E+02) 
2.G857E+02) 
8. 9783E+01 ) 


State vector (unsealed): x T = ( a , 0 T » u f > > £1 . £2 > £1 , (2 ) 
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Table 3.1, concluded 


El GENU ALLIES 


5: 

-.7264 


J 

8.7584 

S: 

-.7264 

- 

J 

8.7584 

7: 

-.4564 

+ 

J 

21.3506 

8: 

-.4564 

- 

J 

21.3506 


EIGENUECTORS: (MAGNITUDE, PHASE (DEGREES) ) 


C l.OOOOE+OOf 
( 7.7404E+00, 
C 6.2283E-03, 
C 8.8075E-01, 
( 8.8783E+00, 
( 2.4390E-01, 
( 7.8027E+01, 
( 2.1435E+00, 


5 

0) 

6.8135E+01) 
7 . 1265E+01 ) 
-2.6606E+01 ) 
-2.3329E+02) 
-2* 1916E+02) 
-1.3855E+02) 
-1.2442E+02) 


( 1.0000E+00. 
( 7.7404E+00, 
( 6.2283E-03, 
( 8.8075E-01, 
( 8.8783E+00, 
( 2.439QE-01. 
( 7.8027E+01, 
( 2# 1435E+00, 


( 1.0000E+00, 
( 2.0127E+01# 
( 2.6446E-Q3, 
( 9. 4247E-01 , 
( 3.3014E+00, 
( B.0949E+01, 
( 7 • 0502E+01 , 
( 1.3016E+03, 


7 

0) 

6.218QE+01) 

7.2886E+01) 

-2.9045E+01) 

2.3713E+02) 

2.8848E+02) 

-3.1647E+01) 

1.9709E+01) 


( 1 . Q0Q0E+00, 
C 2.0127E+01, 
( 2.6446E-03, 
( 9.4247E-01, 
( 3.3014E+00, 
( 6.0949E+01, 
( 7 . 0502E+01 > 
( 1.3016E+03, 


6 

0) 

•6.8135E+01 ) 
7 • 1265E+01 ) 
2.6606E+01) 
2.3329E+02) 
2* 1916E+02) 
1.3855E+02) 
1.2442E+02) 


8 

0) 

•6.218QE+01 ) 
■7.288BE+01) 
2.9045E+01 ) 
■2.3713E+02) 
■2.8848E+02) 
3, 1647E+01 ) 
T . 9709E+01 ) 



MAGNITUDE MAGNITUDE MAGNITUDE MAGNITUDE 
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IMPULSE RESIDUE MRGNITUDES 
(DUE TO PILOT INPUTS) 
CONE 2 OPEN LOOP 
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Figure 3.3 

Impulse Residue Magnitudes - Open-Loop 
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Figure 3.4, continued 
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Figure 3.4, continued 
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the open-loop impulse residue magnitudes, one can see that the rigid-body 
dynamics are dominated by the short period mode, and the first elastic mode 
contributes much more than the second elastic mode to the overall system 
dynamics. Therefore, for this synthesis, choose the short period and first elas- 
tic eigenvalues as the ones to be specified. 

The desired short period complex conjugate eigenvalue pair will be placed 
at the same location as in the “more rigid” baseline configuration. In an 
attempt to reduce the contribution of the first elastic mode to the overall sys- 
tem dynamics, the first elastic complex conjugate eigenvalue pair will be placed 
at a location such that the mode’s natural frequency is the same as the open- 
loop natural frequency and the mode has increased damping. 

Therefore choose: 

X dj8 = X djp = - 1.49 ± j 2.37 = K= 2.799 (rad/s) , *= 0.532) 
x d 34 = x d t] = " 1-758 ± j 8.611 = (w n = 8.789 (rad/s) , f= 0.200) 


Desired Eigenvectors. Specify desired eigenvectors for the short period and first 
elastic modes. Choosing proper desired eigenvectors is very important. For this 
example, they must be chosen such that the closed-loop system will retain an 
“aircraft like” eigenspace. 

In order to make the flexible aircraft behave more like a rigid aircraft, the 
desired eigenvectors are chosen in such a way as to eliminate the elastic mode’s 
contribution to the rigid-body pitch angle and rigid-body pitch rate outputs. 
Therefore, choose the elements of the desired eigenvectors so as to reduce the 
impulse residue magnitudes associated with the flexible modes in the rigid-body 
degrees of freedom, and to reduce the impulse residue magnitudes associated 
with the rigid-body modes in the flexible degrees of freedom. One way this can 
be accomplished is by choosing the elements of the short period eigenvector 
that are associated with the (scaled) states £ 2 > an d £2 1° h ave small mag- 
nitudes. The desired short period eigenvector element associated with d T is 
chosen to be unit magnitude as a reference value. 

The elements of the desired first elastic eigenvector that are associated 
with the (scaled) rigid-body states u f , 0 r , 7, and 0 T are chosen to have small 
magnitudes. The element associated with the state is chosen to be unit mag- 
nitude as a reference value. The phase relationships of these desired 
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eigenvectors are arbitrarily chosen to be similar to those of the “more rigid” 
baseline configuration’s short period and first elastic eigenvectors, respectively. 
Values chosen for the desired eigenvectors are given in Table 3.2. In ttxis table, 
the * denotes elements that are not weighted in the cost function. Therefore, 
the des|red value for these elements is taken as arbitrary. The weighting on 
each of the placed elements is chosen to be unify. 


Calculation of w ; 's. Values for w ; i — 1,...,/ were calculated using Equation 
(3.10) to be: 

w ? X> = wj = ( -0.0550 - j 0.0288 , -1.606 - j 0.7783 ) 

W3 T 4 = w£ = ( 0.0301 + j 0.0174 , 0.0221 - j 1.099 ) 

Achievable Eigenvectors. Tfie achievable short period and first elastic eigenvec- 
tors were calculated usjng Equation (3.7). These eigenvectors are given in 
Table 3.3- Note that in this table the eigenvectors have been scaled such that 
the first element is unity. 

Closed-Loop System Analysis 

The gain matrix yielding the achievable eigenspace in the closed-loop sys- 
tem and the closed-loop system matrices for this control synthesis are given in 
Appendix A. 3. The augmented system eigenvalues, the eigenvalues of A aug are: 

X sp = - 1.49 ± j 2.37 = (w n = 2.799 (rad/s) , f= 0.532) 


\ phl = + 0.0117 
\ ph2 = + 0.3273 


~ 1.758 ± j 8.611 = (w n = 8.789 (rad/s) , f= 0.200) 
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Table 3.2 

Desired Eigenvectors - Example 3.1 


DESIRED EIGENUECTORS 
EIGENUALUES 

i: -1.4000 + J 2.3700 
2: -1.4000 - J 2.3700 
3: -1.7577 + J 8.6109 
4: -1.7577 - J 8.6109 


EIGENUECTORS: (MAGNITUDE* PHASE ( DEGREES ) ) 


( * 

( 1 . 0000E+0Q 

( * 

( * 

( 5. 000QE-03 

( 5.0000E-04 

( 1.0000E-02 

( 5. 00O0E-03 


( 1.0000E-03* 

C 1.0G00E-02. 
( 1.0000E-04, 

( 1.0000E-03, 

( * * 

( * * 

( 1.0GQ0E+00. 

( * * 


1 

* ) 

1 . 0000E+02) 
* ) 

* ) 

3. 0000E+01 ) 
-1.7000E+02) 
1.5000E+02) 
-5. 0G00E+01 ) 


3 

0) 

6 . 00006 + 01 ) 
7 . OOOOE+Ol ) 
-3.0000E+01) 
* ) 

* ) 

1.4000E+02) 
* ) 


( * 

( 1.0000E+00 

C * 

( * 

( 5.0000E-03 

C 5.0000E-04 
( 1.0000E-02 

( 5.0000E-03 


C 1.0000E-03 
( l.QOQOE-Oa 
( 1.0000E-04 

( 1.0000E-03 

( * 

( * 

( l.OQOOE+OO 
( * 


2 

* ) 

-1 • Q000E+02) 
* ) 

* ) 

-3. 0000E+01 ) 
1 . 7000E+02) 
-1.5000E+02) 
5.0000E+01) 


4 

0) 

-6. 0000E+01 ) 
-7.0000E+01) 
3 . 0000E+01 ) 
* ) 

* ) 

-1.4000E+Q2) 
* ) 


State vector (unsealed): x T = ( oc , 0 r , Uf , 0 r , , £ 2 , £i , £2 ) 
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Table 3.3 

Achievable Eigenvectors - Example 3.1 


ACHIEUABLE EIGENVECTORS 
EIGENVALUES 

1: -1.4900 + J 2.3700 
2: -1.4900 - J 2.3700 
3: -1.7577 + J 8.G109 
45 -1.7577 - J 8. 6109 


E I GENUECTORS : C M AGN I TUDE , PHASE ( DEGREES ) ) 


1 2 


( 

l.OOOOE+OOt 

0) 

( 

1.0000E+00, 

0) 

( 

2.5129E+00, 

9.8270E+01 ) 

( 

2, 5129E+00* 

-9.8270E+01 ) 

( 

1.9947E-02, 

4. 545GE+01 ) 

( 

1 . 9947E-02* 

-4. 5456E+01 ) 

( 

8. 9765E-01 » 

-2.3888E+01) 

( 

8.97G5E-01, 

2. 3888E+01 ) 

t 

9.0783E-03, 

2. 7450E+01 ) 

C 

9.0783E-03, 

-2.7450E+01) 

c 

5.210GE-02, 

-1 . 7034E+02) 

C 

5.210GE-02, 

1 . 7034E+02) 

( 

2.5414E-02, 

1.49G1E+02) 

( 

2.5414E-02, 

-1.49G1E+02) 

( 

1 . 4587E-01 » 

-4.8185E+01) 

c 

1.4587E-01, 

4. 8185E+01 ) 

c 

1.0000E+00. 

3 

0) 

c 

1 . 0000E+00* 

4 

0) 

( 

1.9832E+00, 

9.2307E+01) 

c 

1.9832E+00, 

-9.2307E+01) 

( 

3. 8G22E-03* 

7 . G548E+01 ) 

( 

3.8G22E-03, 

-7.6548E+01) 

( 

2.256GE-01, 

-9.2295E+00) 

c 

2. 2566E-01 > 

9. 2295E+00) 

( 

2.3543E+01* 

7 • 01 lOE+Ol ) 

( 

2.3543E+01* 

-7. 01 10E+01 ) 

( 

1.0039E-01, 

3. G417E+01 ) 

c 

1 . 0039E-01 » 

-3.G417E+01) 

c 

2.0G91E+02, 

1 . 71G5E+02) 

( 

2. 0691E+02* 

-1.71G5E+02) 

( 

8.8227E-01* 

1.3795E+02) 

( 

8.8227E-01, 

-1.3795E+02) 


State vector (unsealed): & T = ( a , 6 T , u f , 6 r , , ( 2 > £i > £2 ) 
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X^ 2 = - 0.393 ± j 20.85 = (w n = 20.85 (rad/s) , ?= 0.019) 


The eigenvectors of A aug are given in Table 3.4. The augmented system resi- 
due magnitudes, for u p = impulse, for each output are given in Figure 3.5. 
The augmented system output time responses due to a unit step input in u p , 
with zero initial conditions, are shown in Figure 3.6. 


Comments on Example 3.1 

This example illustrates some of the properties of the DEA synthesis tech- 
nique. As can be seen by examining the desired and achievable eigenvector ele- 
ments. all the desired elements were not exactly obtained (as the theory states). 
But overall, the achievable eigenspace is close to the desired eigenspace. 
Choosing the desired eigenvectors to uncouple the rigid and flexible modes has 
achieved the desired synthesis result. This is reflected in the closed-loop system 
impulse residue magnitudes. The control law has reduced the contributions of 
the first elastic mode in the rigid-body pitch, rigid-body pitch rate and plunge 
acceleration outputs. The impulse residue magnitudes show that the aircraft’s 
flexible modes are not significantly affecting the 0 r and 9 T outputs. Therefore, 
in this respect, this control law is making the augmented aircraft behave more 
like a rigid aircraft. 

One disadvantage that is readily apparent from this example is that this 
technique does not guarantee a stable closed-loop system. Experience has 
shown that the choice of the desired eigenspace can significantly affect the 
movement of the (n - /) unplaced poles. In this example, the desired eigen- 
space has caused the unplaced phugoid mode to split into two unstable real 
eigenvalues (0.0117 and 0.3273). As was shown in the DEA development, with 
this technique only / of the closed-loop system eigenvalues are placed to desired 
locations. The designer has no direct control over the remaining (n - /) eigen- 
values and no way of predicting a priori their directions of movement. 

At this point the designer has two options: 1) increase the number of 
measurements and thereby place a larger number of eigenvalues or 2) use phy- 
sical insight into the problem and conventional control techniques to achieve 
stability by making individual loop gain adjustments and hopefully not disturb 
the achieved eigenspace too much. Option one is the simplest from a theoreti- 
cal viewpoint, but difficult and costly from a practical viewpoint. This method 
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Table 3.4 

Closed- loop Eigenvectors - Example 3.1 


DEA 4 MEASUREMENTS 


EIGENUALUES 

l: -1.4900 + J 2.3700 
2: -1.4900 - J 2.3700 
3S .3273 + J 0 
4S .0117 + J 0 


E I GENUECTORS : (MAGNITUDE, PHASE ( DEGREES ) ) 




i 



2 

c 

1.0000E+00, 

0) 

C 

1.0000E+00* 

0) 

c 

2. 5129E+0Qf 

-2.6I73E+02) 

( 

2.5129E+00, 

2* 6173E+02) 

c 

1.9947E-02, 

“3. 1454E+02) 

c 

i,9947E-02, 

3* 1454E+02) 

( 

9.97G5E-01, 

-2.3888E+01) 

( 

8.9765E-01* 

2, 3888E+01 ) 

( 

9 . 0783E-Q3» 

-3. 3255E+02) 

( 

9.0783E-03, 

3. 3255E+02) 

c 

5.210Gp-02* 

-1.7034E+02) 

c 

5.210GE“02» 

1.7034E+02) 

( 

2.5414E-Q2, 

~2. 1039E+O2) 

( 

2.5414E-02* 

2, 1039E+0H) 

( 

1 .4587E-01 » 

-4* 8135E+0 1 ) 

( 

1 •4587E“01> 

4. 8i85E+01 ) 



3 



4 

( 

1.0000E+00* 

0) 

( 

1 . OOOOE+OO. 

0) 

( 

i . 1652E+Q0* 

0) 

( 

1.3844E-02, 

0) 

c 

4, 1763E-01 * 

i • 8000E+02) 

( 

1.8182E+00. 

-1.8000E+02) 

c 

3* 5600E+QQ* 

0) 

c 

1.1878E+00, 

0) 

( 

6, 7394E-01 » 

0) 

( 

4.1801E-01. 

0) 

i 

4. 4342E-0H* 

1 • 8000E+02) 

( 

5. 1817E-02* 

-1.8000E+02) 

c 

2. 2059E-Q1 » 

0) 

c 

4.872QE~03> 

0) 

( 

1.4514E-02. 

i • 8000E+02) 

( 

6.0393E-04* 

-1.8000E+02) 


State vector (unsealed): x T = ( a , 8 r , u f , 0 r , , £ 2 ) 
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Table 3.4, concluded 


El GENU ALLIES 

5: -1.7577 + J 8.G109 
6S -1.7577 - J 8. BIOS 
7: -.3931 + J 20.8439 
8: -.3931 - J 20.8499 


EIGENUECTORS: ( MAGNITUDE* PHASE ( DEGREES ) ) 


C l.OOOOE+OO* 
C 1.9832E+00, 
C 3. 86222-03, 
C 2.25G6E-01. 
( 2. 3543E+01 * 
C 1 . 0039E-01 > 
C 2. 0691E+02t 
( 8. 8227E-01 » 


( 1.0000E+00. 
( 1.0405E+Q2. 
( 8.6732E-03, 
( 4.9898E+Q0, 
C 9.2935E+01 » 
( 2.4828E+02, 
( 1.9380E+03. 
C 5. 1776E+03* 


5 

0) 

-2.67G9E+02) 
-2.8345E+02) 
-9. 2295E+00 ) 
~2. 8989E+02) 
-3.2358E+-02) 
-1.8835E+02) 
-2.2205E+02) 


7 

0) 

-2.G100E+02) 
-2.1085E+02) 
6. 7918E+01 ) 
-1.4567E+01) 
-5.2830E+01) 
7.6513E+01) 
3.S250E+01) 


( 1.0GGGE+00, 
( 1.3832E+00. 
C 3.8622E-03, 
( 2. 256GE-01 » 
( 2.3543E+01 > 
C 1.0039E-01, 
C 2.0691E+02, 
( 8.8227E-01, 


C 1.0000E+G0* 
C 1.Q405E+02. 
( 3.S732E-03, 
( 4.9893E+00, 
C 9.2935E+01* 
( 2.4828E+02, 
( 1.9380E+03, 
( 5.1776E+03, 


6 

0) 

2. 67S9E+02) 
2.8345E+02) 
9. 2295E+00 ) 
2. 8989E+02) 
3.2358E+02) 
1.8835E+02) 
2 . 2205E+02) 


8 

0) 

2. 0 100E+02) 
2.1035E+02) 
-6.7918E+01) 
1.45G7E+01) 
5.2830E+01) 
-7.6513E+01) 
-3.8250E+01) 



MAGNITUDE MAGNITUDE MAGNITUDE MAGNITUDE 
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IMPULSE RESIDUE MAGNITUDES 


(DUE TO PILOT INPUTS) 
DEA 4 MEASUREMENTS 
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Figure 3.5 

Impulse Residue Magnitudes - Example 3.1 
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DEfl 4 MERS 




Figure 3.6 

Output Time Responses - Example 3.1 
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Figure 3.6, continued 
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is demonstrated in Example 3.2. Option two can be useful to the designer, but 
it requires enough physical insight into the problem to know which loop gains 
to adjust and by how much. Such a situation is not undesirable but, many 
times this requires an iterative control synthesis strategy. However, this option 
does yield a control law that does not require additional sensors. This indivi- 
dual loop gain adjustment technique is illustrated in Example 3.3. 


Example 3.2 : Six Measurements 

This example will illustrate the DEA control synthesis technique when the 
model is eighth order, with two controls, but now has six measurements avail- 
able for feedback. The design objectives for this example are the same as those 
in Example 3.1. Now the designer can place six eigenvalues. This synthesis 
will take advantage of this extra freedom by placing the phugoid poles (which 
were unstable in Example 3.1) to more desirable locations. The measurement 
vector for this example will consist of the four measurements used in Example 
3.1, plus two velocity measurements. The measurement vector is 

z T =(0 t ,0 t , n Zp , n Zt , v Zp , v z< ) 

where 

0 t = total pitch angle (radians) 

0 t — total pitch rate (radians/sec) 

n Zp = plunge acceleration at the cockpit (g' s) 

n Zt = plunge acceleration at an aft fuselage station (g' s) 

v Zp = velocity at the cockpit = f n Zp (n Zp in g's) 

v Zt = velocity at an aft fuselage station = fn Zt (n Zt in g's) 


Specification of Synthesis Parameters 

Desired Eigenvalues. The desired short period and first elastic eigenvalue loca- 
tions are chosen to be the same as in Example 3.1. The phugoid poles will be 
specified to be at the same locations as in the baseline configuration. Therefore 
choose: 
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X d ll2 = X d 5p = “ 1-49 ± j 2.37 = (w n = 2.799 (rad/s) , f= 0.532) 
X d 3 , 4 = X d ph = ~ 0.0015 ± j 0.0672 = (w n = 0.067 (rad/s) , 0.022) 


x d 6 , 6 = x d fi = “ 1-758 ± j 8.611 = (w n = 8.789 (rad/s) , 0.200) 


Desired Eigenvectors. The desired short period and first elastic eigenvectors are 
specified to have the same values as in Example 3.1. The desired phugoid 
eigenvector is chosen using the same strategy as in Example 3.1. The elements 
of this eigenvector that are associated with the (scaled) elastic states are chosen 
to have small magnitudes. The element associated with the the (scaled) state 
tif is choseii to be unit magnitude as a reference value. The phase relationships 
of this desired eigenvector are arbitrarily chosen to be similar to that of the 
“more rigid" baseline configuration’s phugoid eigenvector. The desired eigen- 
vectors for this example are given in Table 3.5. 


Calculation of W;’s. Because the desired short period and first elastic eigen- 
values and eigenvectors are chosen to be the same as in Example 3.1, the 
values for Wj T 2 = w S p and w^ 4 = w^ are also the same. Values for w^ 6 , the W; 
associated with the desired phugoid eigenvector, were calculated to be: 

W 5 T 6 = w p T h = (0.057 + j 8.56 x 10~ 4 , -1.553 + j 6.95 x 10" 3 ) 

Achievable Eigenvectors. The achievable eigenvectors u . ai2 = v.^ and 

are the same as those obtained in Example 3.1. The achievable eigenvectors 
for this example are given in Table 3.6. Again, the eigenvectors have been 
scaled such that their first element is unity. 


Closed-Loop System Analysis 

The gain matrix to obtain the achievable eigenspace in the closed-loop sys- 
tem and the augmented system matrices are given in Appendix A.3. The sys- 
tem closed-loop eigenvalues, the eigenvalues of A aug are: 
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Table 3.5 

Desired Eigenvectors - Example 3.2 


DESIRED EIGEMUECTORS 


EIGEMUALUES 


IS 

-1.4900 

+ 

J 

2.3700 

2: 

“1.4900 

- 

J 

2.3700 

32 

-.0015 

+ 

J 

.0672 

42 

-.0015 

- 

J 

. 0G72 

52 

-1.7577 

+ 

J 

8.6109 

6: 

-1.7577 

- 

J 

8.6103 


E I GENUECTORS : C MAGN I TUDE , PHASE C DEGREES ) ) 


( * > 

( 1 • 0000E+00? 

C * . 

( * » 

( 5. 0000E-03* 

C 5. 0000E-04 » 
( 1.00Q0E-02* 

( 5. 0000E-03* 


( * > 

( * > 

( 1.0000E+00, 

C * 

( 1.0000E-02* 

C 5. 0000E”03f 
C 5.0000E-04, 
( 1.0000E-04, 


C 1.0000E-03, 
C 1.0000E-02. 
C 1.0000E-04, 
( 1.0000E-03, 
( * > 
( * » 
( 1 . OOOOE+OOt 
( * * 


1 

* ) 

1.0000E+02) 
* ) 

* ) 

3. OOOOE+Ol ) 
-1.7000E+02) 
1 . 5000E+02) 
~5. OOOOE+Ol ) 


3 

* ) 

* ) 

1.8000E+02) 
* ) 

1 . 0000E+00) 
1 • 8000E+02) 
3. 1000E+01 ) 
-3. 0000E+01) 


5 

0) 

6. 0000E+01) 
7 . 0000E+01 ) 
-3. OOOOE+Ol) 
* ) 

* ) 

1 . 4000E+02) 
* ) 


C * * 
( l.OOOOE+OO. 
C * t 
( * 

( 5.0000E-03, 
( 5.0000E-04, 
( 1.0000E-02, 
C 5.0000E-03, 


( 1 . 0000E-03* 
( 1.0000E-02, 
C 1.000QE-04, 
( 1 . 0000E-03* 
( * » 
( * » 
C l.OOOOE+OO, 
( * , 


2 

* ) 

-1.G000E+02) 
* ) 

* ) 

-3. OOOOE+Ol ) 
1.7000E+02) 
-1.5000E+02) 
5. OOOOE+Ol) 


G 

0) 

-6. OOOOE+Ol) 
-7. OOOOE+Ol) 
3. OOOOE+Ol) 
* ) 

* ) 

-1 . 4000E+02) 
* ) 


l.OOOOE+OO, 
* , 
1.0000E-02, 
5. 0000E-03, 
5.0000E-04, 
1.0000E-04, 


4 

* ) 

* ) 

-1.8000E+02) 
* ) 

-1 . OOOOE+OO ) 
~1 . 8000E+02) 
-9. 1000E+01 ) 
3, OOOOE+Ol ) 


State vector (unsealed): xJ = ( a , 0 T , u f , 0 r , & , , ft , ft ) 
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Table 3.6 

Achievable Eigenvectors - Example 3.2 


A CHIEUABLE EIGENUECTORS 


EIGENUALUE5 


IS 

-1*4900 

+ 

J 

2*3700 

2: 

-1.4900 

- 

J 

2.3700 

3: 

-.0015 

+ 

J 

. 0G72 

4S 

-.0015 

- 

J 

• 0G72 

5: 

-1.7577 

+ 

J 

8.G109 

G: 

-1.7577 

- 

J 

8.6103 


EIGENUECTORS: (MAGNITUDE* PHASE ( DEGREES ) ) 


C I.OOOOE+OOp 
C 2.5129E+00* 
C 1.9947E-02, 
( 8.9765E-01, 
C 9.0783E-03, 
( 5.2106E-02* 
( 2.5414E-02, 
( 1 . 4587E-01 * 


C 1.0000E+00, 
( 3.2477E-01* 
( 2.3926E+00, 
C 4.8317E+00* 
C 2.3432E-02, 
( 5.1519E-02* 
( 1.5750E-03, 
C 3.4G29E-03* 


( 1.0000E+00, 
( 1 •9832E+00* 
C 3.8622E-03, 
( 2.25G6E-01* 
( 2.3543E+01* 
( 1.0039E-01, 
( 2.0G91E+02* 
( 8.8227E-01* 


i 

0) 

9. 8270E+01 ) 
4. 545GE+01 ) 
-2. 3888E+01 ) 
2.7450E+01) 
-I . 7034E+02) 
1 . 49G1E+02) 
-4. 8185E+01 ) 


3 

0) 

1 . 71G2E+02) 
-1.7913E+02) 
8. 0354E+01 ) 
2, 0432E+00 ) 
1.8023E+02) 
9. 3305E+01 ) 
-8.8511E+01) 


5 

0) 

9. 2307E+01 ) 
7 • G548E+01 ) 
-9.2295E+00) 
7 . 01 10E+01 ) 
3.G417E+01 ) 
1.71B5E+02) 
1.3795E+02) 


( 1.0000E+00* 
C 2.5129E+00* 
( 1.9947E-02* 
( 8.9765E-01* 
C 9.0783E-03, 
C 5.210GE-02* 
( 2. 5414E-02* 
C 1.4587E-01* 


( 1.0000E+00* 

( 3.2477E-01, 

C 2.392GE+00* 
C 4.8317E+G0* 
( 2.3432E-02* 

C 5. 1519E-02* 
( 1.5750E-03* 

C 3.4G29E-03* 


2 

0) 

-9.8270E+Q1) 
-4.5456E+01) 
2. 3888E+01 ) 
-2.7450E+01) 
1.7034E+02) 
-1.49G1E+02) 
4. 8185E+01 ) 


4 

0) 

-1.71G2E+02) 
1.7913E+02) 
-8. 0354E+01 ) 
“2. 0432E+00) 
-1.8023E+02) 
-9.3305E+01) 
8*851 1E+01 ) 


( 1.0000E+00f 
( 1.9832E+00* 
C 3.8622E-03* 
C 2. 2566E-01* 
( 2. 3543E+01 * 
( 1.0039E-01* 
C 2. 0G91E+02* 
C 8.8227E~01t 


6 

0 ) 

-9. 23Q7E+01 ) 
-7.6548E+01) 
9. 2295E+Q0 ) 
-7.0110E+01) 
-3. G417E+01 ) 
-1.71G5E+02) 
-1.3795E+02) 


State vector (unsealed): * T = ( a , 6 T , u f , 6 r , & , & , Zi , ) 
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X sp = - 1.49 ± j 2.37 = (w n = 2.799 (rad/s) , 0.532) 

X ph = - 0.0015 ± j 0.0672 = (w n = 0.067 (rad/s) , f= 0.022) 
X 6 = - 1.758 ± j 8.611 = (w n = 8.789 (rad/s) , f= 0.200) 


X & = - 0.3531 ± j 20.97 = (w D = 20.97 (rad/s) , 0.017) 


The eigenvectors of A aug are given in Table 3.7. The augmented system 
impulse residue magnitudes (expressed as a percentage of the total response 
for each output are shown in Figure 3.7, and the output time responses di” ' 
a unit step in u. p , with zero initial conditions, are given in Figure 3.8. 


Comments on Example 3.2 

For this example, using DEA with the freedom to place six eigenvalues has 
yielded a stable closed-loop system. The short period, phugoid, and first elastic 
eigenvalues have been placed to desired locations. The second elastic eigen- 
values are at an acceptable location. They have remained close to their open- 
loop position. 

The achieved eigenspace is close to the desired eigenspace. The impulse 
residue magnitudes show that the control synthesis has achieved the desired 
result; reducing the residue magnitudes of the first and second elastic modes in 
the rigid-body pitch angle and pitch rate outputs. In this respect, the control 
law is making the flexible aircraft behave more like a rigid aircraft. 

These improved results are obtained at the “cost” of more sensors. This 
synthesis requires six measurements and is therefore more complex and costly 
than the four measurement synthesis of Example 3.1. 

The next example will demonstrate a technique that will stabilize the sys- 
tem without increasing the number of measurements. This technique shall be 
referred to as Individual Loop Gain Adjustment (ILGA). 
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Table 3.7 

Closed-loop Eigenvectors - Example 3.2 


DEA 6 MEASUREMENTS 
EIGENUALUES 

1 : -1.4900 + J 2.3700 
2: -1.4900 - J 2.3700 
3: -.0015 + J . 0672 
4S -.0015 - J • 0672 


E I GEINUECTORS J ( M AGM I TUDE , PHASE ( DEGREES ) ) 


( 1.0000E+00, 
C 2.5129E+00, 
( 1.9947E-02, 
( 8.97G5E-01, 
( 9.0783E-03, 
( 5.210GE-02, 
C 2.5414E-02, 
C 1.4587E-01, 


C 1 • OOOOE+OOf 
C 3.2477E-01, 
C 2. 3926E+00* 
( 4.8317E+00, 

( 2.3432E-02, 

C 5.1519E-02* 
( 1.5750E-03, 

C 3.4629E-03, 


1 

0) 

9.8270E+01) 
4.5456E+01) 
-2.3888E+01) 
2. 7450E+01) 
1 • 8966E+02) 
1 . 4961E+02) 
-4.8185E+01) 


3 

0) 

1.7162E+02) 
1 . 8087E+02) 
8.0354E+01) 
2. 0432E+00) 
1 . 8023E+02) 
9.3305E+01) 
-8.8511E+01) 


( 1. Q000E+00* 
C 2.5129E+00* 
( 1.9947E-02, 
( 8.97G5E-01* 
( S.0783E-03, 
C 5.2106E-02* 
C 2.5414E-02* 
C 1 .4587E-01 » 


2 

0) 

-9.8270E+01) 
-4.545GE+01 ) 
2.3888E+01 ) 
-2.7450E+01) 
-1.89p6E+02) 
-1.49G1E+02) 
4.8185E+01 ) 


( 1 * 0Q00E+00* 

( 3.2477E-01, 

( 2. 3926E+00* 

( 4.8317E+00, 

( 2.3432E-02, 

C 5.1519E-02, 
( 1.5750E-03, 

( 3.4G29E-03, 


4 

0) 

-1.71G2E+02) 
-1.8087E+02) 
-8.0354E+01) 
-2.0432E+00) 
-1.8023E+02) 
-9.3305E+01) 
8.851 1E+01 ) 


State vector (unsealed): x T = ( a , 0 T , u f , B r , ( x , £> , & , ) 
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Table 3.7, concluded 


EIGENVALUES 


5: 

-1.7577 


J 

8.G108 

6: 

-1.7577 

- 

J 

8.6109 

7: 

-.3531 

+ 

J 

20.9739 

8: 

-.3531 

- 

J 

20.9739 


EIGENUECTORS: (MAGNITUDE* PHASEC DEGREES ) ) 


C I . 0000E+00* 
( 1.9832E+00, 
( 3.8622E-03, 
( 2.256GE-01* 
( 2. 3543E+01 * 
( 1.0G39E-01* 
( 2. 0691E+02* 
( 8.8227E-01* 


5 

0 ) 

9. 2307E+Q1 ) 
7.6548E+01) 
-9.2295E+00) 
7.0110E+01 ) 
3.6417E+01 ) 
1.7165E+02) 
1.3795E+02) 


C 1.0000E+Q0. 
( i . 8194E+0 1 9 
( 2.6346E-03* 
( 8.G732E-01, 
( 5.1347E+Q0* 
( 2.58G3E+01* 
( 1 . 0771E+02* 
( 5.4252E+02* 


7 

0 ) 

1.0645E+02) 
9, 7280E+0 1 ) 
1.5430E+01) 
1.7445E+01) 
2. 1957E+02) 
1.0841E+Q2) 
-4. 3466E+0 I ) 


( 1 . 0000E+00* 
C 1 • 9832E+00* 
( 3.S622E-03* 
( 2.2566E-01, 
( 2. 3543E+Q 1 * 
( 1.0039E-01, 
( 2.0G91E+02, 
( 8.8227E-Q1* 


C l.OOOOE+OO* 
( 1.8194E+Q1. 
( 2.6346E-03* 
( 8.S732E-01 9 
( 5. 1347E+00* 
( 2.58S3E+01. 
( 1 . 0771E+02* 
( 5.4258E+02* 


6 

0) 

9* 2307E+01 ) 
•7.G548E+01) 
9.2295E+00) 
7 . 01 10E+01 ) 
3. 6417E+01 ) 
■1.71G5E+02) 
1 . 3795E+02) 


8 

0 ) 

1 . 0645E+02) 
■9. 7280E+C1 ) 
1.5490E+01) 
1 . 7445E+01 ) 
2* 1957E+02) 
1 . 0841E+02) 
4.3466E+01) 



MAGNITUDE MAGNITUDE MAGNITUDE MAGNITUDE 
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IMPULSE RESIDUE MRGNITUDES 
(DUE TO PILOT INPUTS) 
DER 6 MEASUREMENTS 
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Figure 3.7 

Impulse Residue Magnitudes - Example 3.2 
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Figure 3.8 

Output Time Responses - Example 3.2 
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Figure 3.8, continued 
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Example 3.3 : Individual Loop Gain Adjustment 

This example will illustrate the use of an individual loop gain adjustment 
technique to stabilize the unstable system that resulted in Example 3.1. This 
method relys upon the designer’s physical insight into the problem. The 
designer must choose which loop gains to vary to stabilize the system. If only 
a small change in loop gain is needed to stabilize the system, hopefully the 
resulting eigenspace will not be disturbed too much from the desired eigen- 
space. 

This method involves choosing which loop gain to vary based upon the 
transfer function characteristics for this measurement to input and then using 
the Root Locus method to determine the change in the loop gain required to 
obtain acceptable eigenvalue locations. The new closed-loop eigenspace must 
then be checked to see how close it is to the original achievable eigenspace. 

This example will start with the results from Example 3.1. The closed- 
loop system eigenvalues obtained were: 

X sp = - 1.49 ± j 2.37 = (w n = 2.799 (rad/s) , f= 0.532) 


X phl = + 0.0117 
X ph2 = + 0.3273 

X ?i = - 1.758 ± j 8.611 = (w n = 8.789 (rad/s) , f= 0.200) 


X ?2 = - 0.393 ± j 20.85 = (w n = 20.85 (rad/s) , f- 0.019) 

As can be seen, it is the phugoid mode that is unstable. This mode “shape” is 
predominantly pitch attitude, 6 and forward velocity, u f , and therefore the pole 
locations should be affected more by the elevator than the control vane. There- 
fore, as a first attempt, the effect of varying the <5 e -to-0 t loop gain will be 
explored. The first step of this process is to calculate the 0 t -to-<5 e transfer func- 
tion. The zero’s of the 0 t -to-<5 e transfer function were calculated using the 
method of Sandburg and So [14] to be: 



58 


z j = -0.00917 z 2 = -0.745 
z 3,4 = ~0-856 ± j 5.365 
z 5 6 = -0.864 ± j 29.8 


Therefore, the 0 t -to-<5 e transfer function is 

A = k(s - Zj)(s - z 2 )(s - z 8<4 )(s ~ z 8iB ) 

( s “ \hl)( s " ^ph2)( s ~ ^sp)( s “ k^)(s ~ X^) 

Using the Root Locus method, the 0 t -to-<5 e feedback gain was varied in a-?, 
attempt to stabilize the system (see Figure 3.9 and 3.10). At a gain of, AG$ $ 
= 0.28 (dimensionless) the system eigenvalues are: 

X sp = - 1.243 ± j 2.907 = (w n = 3.16 (rad/s) , f= 0.393) 


X ph = - 0.0253 dfc j 0.068 = (w n = 0.073 (rad/s) , ?= 0.345) 


X ?i = - 1.844 ± j 9.100 = (w n = 9.28 (rad/s) , f= 0.199) 


X ?2 = - 0.3589 ± j 20.69 = (cu n = 20.70 (rad/s) , f= 0.017) 


The closed-loop system gain to achieve this eigenvalue placement is 

G' = G + 




o 

!o‘ 

o' 

;o; 

r o" 

io. 


where G is the gain matrix obtained in Example 3.1 

This modified gain matrix and the augmented system matrices are given in 
Appendix A. 3. The eigenvectors of the augmented system are given in Table 
3.8. The augmented system impulse residue magnitudes (normalized for plot- 
ting) for each output are shown in Figure 3.11, and the output time responses 
due to a unit step in u p , with initial conditions x(0) = 0, are shown in Figure 
3.12. 
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□ Closed-loop Pole 
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Figure 3.9 ▼ 

0 t /<5 e Transfer Function Root Locus 
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O Transfer Function Zero 

X Transfer Function Pole 

□ Closed-loop Pole 



Figure 3.10 

Root Locus, Detail near Origin 
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Table 3.8 

Closed-loop Eigenvectors - Example 3.3 


DEA ILGA 
EIGENUALUES 

l: -1.2432 + J 2.9073 
2: -1.2432 - J 2.9073 
3: -.0253 + J . 0688 
4: -.0253 - J .0G88 


EIGENUECT0R5: (MAGNITUDE, PHASE (DEGREES) ) 


1 


( 

1 . OOOOE+OO, 

0) 

( 

1.0000E+00* 


( 

3. G365E+00 , 

9. 1774E+01 ) 

( 

3. 03G5E+00, 

-9 

( 

1.8373E-02, 

5 . 5431E+01 ) 

( 

1.8373E-02, 

-5 

( 

9.6030E-01. 

-2.1378E+01) 

( 

9.G030E-01, 

2 

( 

1.3229E-01, 

-6. 2820E+00) 

( 

1.3229E-01, 

G 

( 

G.1202E-02, 

1 . 86G2E+02) 

( 

G.1202E-02, 

-1 

( 

4.1830E-01. 

1 . 0G87E+02) 

( 

4.1830E-01. 

-1 

( 

1.9352E-01, 

-6. 0227E+01 ) 

( 

1.9352E-01, 

G 

( 

K0000E+00. 

3 

0) 

( 

1.0000E+00. 


( 

3. G050E-01 , 

-1.4212E+02) 

( 

3.B050E-01. 

1 

( 

2. 3365E+00 , 

1.8849E+02) 

( 

2. 3365E+00, 

-1 

( 

4. 9170E+00, 

1 . 0771E+02) 

( 

4. 9170E+00, 

-1 

( 

1.3442E+00. 

9. 7165E+01 ) 

( 

1 . 3442E+00, 

-9 

( 

5. 1354E-02, 

-1.6234E+02) 

( 

5.1354E-02, 

1 

( 

9. 8555E-02, 

-1.52GGE+02) 

( 

9.8555E-02, 

1 

( 

3.7651E-03, 

-5. 2157E+01 ) 

( 

3.7G51E-03, 

5 


State vector (unsealed): x. T = ( a , 0 T , u r , 8 X , , £ 2 , , £ 


2 

0) 

1774E+01) 
5431E+01) 
1378E+01 ) 
2820E+G0 ) 
8G62E+02) 
0G87E+02) 
0227E+01) 


4 

0 ) 

4212E+02) 
8849E+02 ) 
0771E+02) 
7165E+01) 
E234E+02) 
52GGE+02) 
2157E+01 ) 


) 
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Table 3.8, concluded 


EIGENUfiLUES 
5: -1.8439 + J 

6: -i . 8439 - J 

7: -.3589 + J 

8: -.3589 - J 


9.0998 

9.0998 

20.6938 

20.6938 


EIGENUECTORS: (MAGNITUDE, PHASE ( DEGREES ) ) 


( 1.0000E+00, 
( 8.7S51E+0G, 
C 5.6834E-03, 
( 9.4296E-01, 
( 1.S359E+01, 
( 4. 0214E-01 > 
( 1 . 7974E+02, 
( 3.7338E+00, 


( 1 • QOOOE+Q0, 
( G.7726E+01, 
( 6.1G38E-03, 
( 3.2723E+00, 
( 4. 9552E+01 » 
( 1 . 2272E+02, 
( 1 • 0256E+03* 
( 2.5398E+03, 


5 

0) 

5.0591E+01) 
5. 0627E+01 ) 
-5.0864E+01) 
1.3417E+02) 
1.3425E+02) 
-1.2438E+02) 
-1.2430E+02) 


7 

0) 

3. 4772E+01 ) 
4. 2737E+01 ) 
3. 0378E+02) 
1 . 9483E+02 ) 
1 . 5740E+02) 
2. 8582E+02) 
2.4840E+02) 


( 1.0000E+00, 
( 8.7551E+00, 
C 5.6834E-03, 
C 9.4296E-01, 
C i • 9359E+01, 
( 4.0H14E-01, 
( i • 7974E+02, 
( 3. 7338E+00, 


( 1 . OOOOE+OO, 
( 6.7726E+01* 
( 6. 1638E-03, 
( 3.2723E+00, 
( 4.9552E+01, 
( 1 . 2272E+02, 
( 1 . 0256E+03, 
( 2.5398E+03, 


6 

0) 

-5. 0591E+01 ) 
-5. 0627E+01 ) 
5. 0864E+01 ) 
-1.3417E+02) 
-1.3425E+02) 
1 . 2438E+02 ) 
1 . 2430E+02) 


8 

0) 

-3.4772E+0I) 
-4.2737E+01) 
-3. 0378E+02) 
-1.9483E+02) 
-1.5740E+02) 
-2.8582E+02) 
-2.4840E+02) 



MAGNITUDE MAGNITUDE MAGNITUDE MAGNITUDE 
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IMPULSE RESIDUE MAGNITUDES 
(DUE TO PILOT INPUTS) 
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Figure 3.11 

Impulse Residue Magnitudes - Example 3.3 
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Figure 3. 12, concluded 
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Comments on Example 3.3 

By using an individual loop gain adjustment technique on the synthesis 
result of Example 3.1, one can stabilize the system without increasing the 
number of measurements. Adjusting the 0 t -to-6 e gain has stabilized the 
phugoid eigenvalues and left the short period, first elastic, and second elastic 
poles relatively close to their original (augmented) locations. 

By examining the closed-loop system eigenvectors one can see that the 
closed-loop system eigenspace obtained using this synthesis has drifted a little, 
but is still close to the achievable eigenspace (that obtained in Example 3.1). 
This eigenspace adjustment is reflected in the closed-loop impulse residue mag- 
nitudes obtained from this synthesis. The synthesis of Example 3.1 yielded a 
control law that had practically eliminated the contribution of the first and 
second elastic modes to the 0 r and 0 T outputs. By comparing the impulse resi- 
due magnitudes from Example 3.1 (Figure 3.5) and this example (Figure 3.11), 
one can see that the eigenspace adjustment has caused a slight increase in the 
contribution of the first elastic mode to the 0 T and 0 t outputs. This can also be 
noted in the 0 t time response. The slight ripple on this response is due to the 
first elastic mode. 

As can be seen from the 0 T and 0 T impulse residue magnitudes, the elastic 
mode's contribution is still small compared to the rigid-body modes contribu- 
tion. Therefore, this synthesis still tends to achieve the desired synthesis objec- 
tives; improving the rigid-body dynamics and making the flexible aircraft 
respond more like a rigid aircraft. 


DEA Conclusions 

DEA is a control synthesis technique for directly determining measurement 
feedback control gains that will yield an achievable closed-loop eigenspace. For 
an observable, controllable system that has n states, m controls, and / measure- 
ments one can determine a gain matrix that will place l eigenvalues to desired 
locations and their l associated eigenvectors as close as possible in a least 
squares sense to desired eigenvectors (assuming / > m). 

Using DEA the designer chooses l desired eigenvalues and their associated 
eigenvectors to achieve a desired synthesis objective (such as adjusting residue 
magnitudes, uncoupling system modes, reducing a certain eigenvalue’s contri- 
bution to an output, response shaping, etc.) The desired eigenspace is used to 
determine the system’s achievable eigenspace. By comparing the desired and 
achievable eigenspace, the designer can obtain an estimate of how close the 
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synthesis will come to the desired synthesis objective. This achievable eigen- 
space is then used to calculate control gains. As the examples demonstrated, 
DEA has the potential to be a useful control synthesis technique. This tech- 
nique can be used alone in an iterative fashion as demonstrated in Examples 
3.1 and 3.2 or, it can be used in conjunction with other control synthesis tech- 
niques (as demonstrated in Example 3.3) to achieve desired results. When 
using DEA in an iterative fashion, the designer synthesizes a control law using 
a small number of measurements. If the synthesis results do not meet the 
desired synthesis objective, the number of measurements is increased to obtain 
more design freedom and a new control law is synthesized. The number of 
measurements is increased until an acceptable synthesis is obtained. 

The examples also presented the two major disadvantages of this tech- 
nique. First, this technique does not always yield improved performance. In 
some cases, it can yield an unstable closed-loop system. Using DEA the 
designer can only place l of the closed-loop poles. The designer has no infor- 
mation concerning the closed-loop locations of the remaining (n - /) eigenvalues 
and no way of predicting beforehand their directions of movement. Experience 
has shown that the choice of the / desired eigenvectors can significantly effect 
the movement of the unplaced poles. 

If all of the states are available for feedback, one could use DEA with 
full-state feedback and place all of the poles of the closed-loop system. An 
example of this using a model that includes actuator dynamics is presented in 
Reference [13]. In this example, by placing all the poles of the closed-loop sys- 
tem, it was possible to obtain acceptable robustness, reasonable actuator 
bandwidths, and achieve the desired synthesis results. This model with actua- 
tor dynamics and the issue of required actuator bandwidths are considered 
further in Chapter IV. 

It is possible for this technique to yield unacceptably high control gains. 
There is no systematic way of trading control energy for eigenspace assignment 
goals. In some designs, exactly obtaining an achievable eigenspace may require 
a large control effort when just getting close would yield satisfactory results. A 
control system that utilizes large control effort may be very susceptible to sen- 
sor noise. 

The next chapter presents a possible alternative to the direct eigenspace 
assignment technique. 
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CHAPTER IV 

LINEAR QUADRATIC EIGENSPACE ASSIGNMENT 


As was shown in Chapter El, Direct Eigenspace Assignment is a very use- 
ful control synthesis technique but it has two major disadvantages. First, there 
is no systematic way to trade eigenspace assignment goals versus gain magni- 
tudes. Second, when full state feedback is not used the (n - /) unspecified sys- 
tem poles can move to unpredictable locations, sometimes even destabilizing 
the system. 

A possible alternative is a synthesis method based upon the work of C.A. 
Harvey and G. Stein [9,10]. This procedure takes advantage of the asymptotic 
modal characteristics of multi-variable linear quadratic regulators. It is based 
upon Harvey and Stein’s procedure for selecting quadratic weighting matrices 
that will asymptotically yield an achievable eigenspace in a linear quadratic 
control synthesis, as the cost function weighting on the control tends toward 
zero. This method shall be referred to as Linear Quadratic Eigenspace Assign- 
ment (LQEA). 


Asymptotic Modal Properties 

In their paper, “Quadratic Weights for Asymptotic Regulator Properties”, 
Harvey and Stein [9] present the following theorem summarizing the asymp- 
totic modal properties of an LQ regulator as a function of the weights Q and 

R. 


Theorem 4.1 


Given the system 


X. = Ax. + Bu 


( system dynamics ) 


(4.1a) 
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L — Hx. ( system responses ) (4-lb) 

where x e R n , u e R m , and r e R m . 

Consider the linear quadratic regulator 

11 = Gx (4.2) 

G = R _1 B T P/p (4.3) 

where P is a solution to the steady state Ricatti equation 

0 = PA + A t P + Q - PBR _1 B T P/p (4.4) 

with cost criterion 

CO 

J = / (r T r + pu T Ru)dt (4.5) 

o 

CO 

= / (x. t Q x + pu T Ru)dt 
o 

where Q = H T H. Assume further that: 

1) Rank (HB) = m 

2) the zero’s of det (H(sl n — A^B) are distinct, have negative real parts, and 
do not belong to the spectrum of A. 

Then, the optimally controlled system has the following properties: 
Asymptotically Finite Modes 

As p tends toward zero there are (n - m) closed-loop system eigenvalues of 
the form 

Hp) -* V i = 1,2, ...,(n - m) (4.6) 


X L r | < oo 
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with associated eigenvectors 

HP) “*■ (^X ~ A) _1 Bwf (4.7) 

where \{ and wf are defined by 

H(Xi f I n - A) _1 Bw/ = 0 (4.8) 


Asymptotically Infinite Modes 

As p tends toward zero there are also m eigenvalues of the form 


Hp) - «j°7 p' 11 

(' 

1 s i°°l < °°l j = 1,2,... ,m 


with associated eigenvectors 

E]{p) -+ Bwj 00 


where Sj°°(p) and Wj°° are defined by 

R = N _t S _2 N _1 

(4.11) 

y t y = (HoB^n-tn-^^b)- 1 

with 

N = [wf° | w 2 °° | • • • | w£>] 

(4.12) 


S = diag [sf°, s$°, ..., s“] 

H = YH 0 


This concludes the Theorem. 
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Discussion of Theorem 

As p tends toward zero (n - m) closed-loop system eigenvalues tend toward 
the placed values 

Xi(/>) \ f i=l,...,(n-m) (4.6) 

These desired finite eigenvalues, are the (n - m) transmission zeros of the 
system given by Equations (4.1a) and (4.1b) [15]. Their associated eigenvectors 
tend toward 

HP) -+ (\ f In _ A) _ 1 Bw / = £ a . (4.7) 

By comparing this result with Equation (3.7) of Chapter III, one can see that 
this equation is just an expression for the achievable eigenvectors of the sys- 
tem. 

The remaining m closed-loop system eigenvalues tend toward infinity in m 
first order Butterworth patterns at a rate of 

XjW - »,'%>''* = Xj°° j = (4.0) 

and their associated eigenvectors tend towards the vectors 

i'j(p) - Bwj 00 = 14 00 (4.10) 

as p tends towards zero. In this development Xj 00 shall be referred to as a 
desired infinite eigenvalue and v?° shall be referred to as a desired infinite 
eigenvector. 

Once values for \(, w/, Sj 00 , and Wj°° are chosen, the Q and R matrices can 
then be determined. The control weighting matrix, R, can be determined 
easily from the equation 

R = N -t S~ 2 N -1 (4.11) 


S = diagfsj 00 ,s 2 °° , • • • s-] 

N = [ wf° | V£ 2 ° | • • • | w“ ] 


where 
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The state weighting matrix, Q, is calculated by first determining achiev- 
able closed- loop system eigenvectors for each \J and w ; f , using Equation (4.7). 
The achievable eigenvectors can be calculated for a given \ ; f and desired eigen- 
vector using the method of Chapter ttt. Concatenate these (n - m) achievable 
eigenvectors columnwise to form 

v = K' 1 14 / 1 ••• | 

Using this matrix, calculate a system projection matrix P n 

P n = I n - ViyTy]" 1 V T (4.13) 


partitioned in the following way 


P„ = 


1*11 1*12 

P i2 P 22 


where P 2 2 ls a nixm submatrix. 


Calculate 


H 0 = [Pa 1 P , t 2 U 


(4.14) 


and 


Y = (H 0 BN)-' 


(4.15) 


The state weighting matrix, Q, can be determined by using 

Q = H t H = (YH 0 ) t (YH 0 ) 




(4.16) 


Closed-loop System 

Using this method, the augmented (closed-loop) system may be written 


i = (A + BG)k + Blip 


(4.17) 


The closed-loop input to state transfer functions are given by 
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2L(s) = [sl n - (A + BG)]- 1 Bu p (s) (4.18) 

The system’s physical output responses are given by 

% = Cx (4.19) 

Therefore, the response transfer functions for the pilot input, u p , are 

X(s) = C (sl„ - (A + BG)] -I Bjip(s) (4.20) 

By selecting the desired eigenspace, one could say that the designer is placing 
the poles and zeros of these transfer functions to desirable locations. This co 
trol synthesis technique is summarized in the following section. 


LQEA Control Synthesis Technique Summary 

The control synthesis is a two part procedure: 

1) Select values for X; f , Sj 00 , and Wj 00 . Use Harvey and Stein’s method to 
determine quadratic weighting matrices Q and R that will asymptotically 
yield the desired closed-loop eigenspace. 

2) Perform an LQ synthesis using various values of control weighting ( p ). 
Select the synthesis that will yield a closed-loop system as close as needed 
to the desired eigenspace. 

Two examples will be presented in the following section to illustrate the 
use of this technique. 


LQEA Examples 

Two examples will be presented to demonstrate some properties of this 
synthesis technique. The two examples differ only in the desired eigenvalue 
and eigenvector for the second elastic mode. The results are intended to 
demonstrate how significantly the desired eigenspace can affect the synthesis 
results. 

The large flexible aircraft model used in Chapter HI assumed “infinitely 
fast” actuators. Both of these examples will also consider this aircraft model, 
but will include elevator and control vane actuator dynamics. These dynamics 
will be modeled as: 
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1 = (I/O 

\ S + (1/t) 


where 

6 [ = actual control surface deflection (radians) 

<5; c = commanded control surface deflection (radians) 

For the control vane actuator, (1.0/r) = 10.0 sec -1 ; for the elevator actuator 
(1.0/r) = 9.0 sec -1 . 


System 

The (unsealed) model that will be used is then of the form: 
x = Ax + Bu (system dynamics) 

y = Cx (system responses) 


11 - Me + Up (total control input) 


He — Gx (feedback control law) 


u p — pilot' s input 


States 


— ( ® ) £l > £2 > £l > £2 > ) ^cv ) 


where 
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a = angle of attack (radians) 

0 T = rigid body pitch rate (radians/sec) 
u f = forward velocity (ft/sec) 

0 T = rigid body pitch angle (radians) 

— mode one generalized deflection (dimensionless) 

£2 = mode two generalized deflection (dimensionless) 

= mode one generalized deflection rate (l/sec) 

= mode two generalized deflection rate (l/sec) 

6 e = elevator deflection (radians) 

8 CV = forward control vane deflection (radians) 

Controls 

a? = ( 6 ,,, ) 

lip" = («.., 0.0 ) 

where 

<5 6c = commanded elevator deflection (radians) 

£ CVc = commanded forward control vane deflection (radians) 


Outputs 

£ T = ( u f » e r > K » 7 , *t » K t n z p ) 


where 
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u f = forward velocity (ft/sec) 

0 t = rigid body pitch angle (radians) 

0 T = rigid body pitch rate (radians/sec) 

7 = flight path angle (radians) 

0 t = total pitch angle (radians) 

0 t = total pitch rate (radians/sec) 

n Zp = plunge acceleration at the cockpit (g 1 s) 

The modified open-loop system matrices can be found in Appendix A.4. 


Example 4.1 

The design objectives for this example are the same as those of Example 
3.1; reducing the contribution of the elastic modes to the rigid-body responses, 
specifically 0 t . For this example dim(x)=n=10 and dim(u)=m=2, therefore one 
may specify eight desired finite eigenvalues/eigenvectors and two desired 
infinite eigenvalues/eigenvectors. 


Specification of Synthesis Parameters 

Desired Finite Eigenvalues. With this technique one has the freedom to specify 
the short period, phugoid, first elastic, and second elastic eigenvalues. Values 
for the locations of the desired finite short period, phugoid, and first elastic 
eigenvalues are chosen to be the same as in Example 3.2 . The second elastic 
eigenvalue is chosen to be a value close to its open-loop location. The damping 
is increased slightly, the natural frequency is taken to be the same as the 
open-loop. 

Xj2 = X s f p = - 1.49 ± j 2.37 = (w n —2.799 (rad/s) , f =0.532) 
x 3 f ,4 = X p f h = “ 0.0015 ± j 0.0672 = (w n =0.067 (rad/s) , $ =0.022) 
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X{ 6 = X£ = - 1.758 ± j 8.611 = K =8.789 (rad/s) , * =0.200) 
X( 8 = X| g = - 0.534 ± j 21.348 = (w n =21.35 (rad/s) , f =0.025) 


Desired Finite Eigenvectors. The desired eigenvectors for the short period, 
phugoid, and first elastic modes are chosen to be the same as specified in 
Example 3.2 . The second elastic eigenvector is chosen to be the same as the 
open-loop second elastic eigenvector. Values chosen for these desired eigenvec- 
tors can be found in Table 4.1. In this table, the * denotes elements whose 
placement is taken as arbitrary. The weighting on each of the specified ele- 
ments is chosen to be unity. 


Achievable Eigenvectors. The achievable eigenvectors are calculated in the 
same way as for the DEA technique. First Equation (3.10) is used to calculate 
values of w ; f for i=l,...,(n-m). Then Equation (3.7) is used to calculate the 
achievable eigenvectors. These achievable eigenvectors are given in Table 4.2. 
In this table, the eigenvectors have been scaled such that the first element is 
unity. 


Desired Infinite Eigenvalues. Values for Sj 00 , which define the “infinite modes”, 
are chosen to maintain a two-to-one actuator bandwidth ratio. Because the 
control vane is physically smaller than the elevator and is expected to play a 
large role in the control of the elastic modes, it is chosen to have the larger 
bandwidth. 

sf = s e f = 0.5 

s 2 °° = s cv = 1-0 


Desired Infinite Eigenvectors. Values for Wj 00 are chosen to achieve decoupled 
actuator eigenvectors. 
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Table 4.1 

Desired Eigenvectors - Example 4.1 


DESIRED EIGENUECTORS 
El GENU ALLIES 

1: -1.4900 + J 2.3700 
2: -1.4900 - J 2.3700 
3: -.0015 + J .0672 
4: -.0015 - J .0672 


EIGENUECTORS: (MAGNITUDE, PHASE (DEGREES) ) 


( * , 

( 1 . OOOOE+OO, 

( * , 

( * , 

( 5. 0000E-03, 

( 5. Q000E-04, 

( 1.0000E-02, 

( 5. 0000E-03, 

( * , 

( * , 


( * 

( * 

( 1. 0000E+00, 

( * , 

( 1.0000E-02, 

( 5.0000E-03, 

( 5.0000E-04, 

( 1.0000E-04, 

( * 

( * 


1 

* ) 

1.0000E+02) 
* ) 

* ) 

3. OOOOE+Ol ) 
-1 . 7000E+02 ) 
1.5000E+02) 
-5. OOOOE+Ol ) 
* ) 

* ) 


3 

* ) 

* ) 

1 • 8000E+02) 
* ) 

1. 0000E+00) 
1.8000E+02) 
9. 1000E+01 ) 
-9. 0000E+01) 
* ) 

• ) 


( * , 

( 1. 0000E+00, 
( * , 
( * , 
( 5. 0000E-03, 
( 5. 0000E-04, 
( 1 . 0000E-02, 
( 5. 0000E-03, 
( * , 
( * , 


( * 

( * 

( 1. 0000E+00 

( * 

( 1.0000E-02 

( 5. 0000E-03 

( 5.0000E-04 

( 1.G000E-04 

( * 

( * 


2 

* ) 

-I • 000OE+02) 
* ) 

* ) 

-3. OOOOE+Ol) 
1 • 7000E+02) 
-1 • 5000E+02) 
5. OOOOE+Ol) 
* ) 

* ) 


4 

* ) 

* ) 

-1.8000E+02) 
* ) 

-1. OOOOE+OO) 
-1 . 8000E+02) 
-9. 1000E+01 ) 
9. OOOOE+Ol) 
* ) 

* ) 


State vector (unsealed): & T = ( a , 0 T , Uf , O r , , £ 2 > £i » £2 > 4 , M 
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Table 4.1, concluded 


EIGENUALUES 

5S -1.7577 + J 8.G109 
G2 -1.7577 - J 8.G109 
7: -.5339 + J 21.3480 
8: -.5339 - J 21.3480 


E I GENUECTORS : C MAGN I TUBE , PHASE C DEGREES ) ) 


C 

( 

( 

( 

( 

( 

C 

c 

c 

( 


1.0000E-03, 
1 • 0Q00E-02* 
1.0000E-04, 
1 . 0000E-03* 
* * 

* , 

1.0000E+00. 


* 


5 

0) 

G.0000E+01) 
7. 0000E+01 ) 
-3. 0000E+01 ) 
* ) 

* ) 

1.4000E+02) 
* ) 

* ) 

* ) 


( 1.0000E-03, 
C 1.0000E-02, 
C 1.0000E-04, 
C 1 .0000E-03# 
( * * 
( * 9 

( 1.0000E+Q0, 
C * » 
( * , 
C * * 


G 

0) 

-6. 0000E+01 ) 
-7. 0Q00E+0 1 ) 
3. 0000E+01 ) 
* ) 

* ) 

-1.4000E+02) 
* ) 

* ) 

* ) 


C 1.000GE+Q0* 
C 2. 0194E+01# 
( 2.G518E-03, 
C 9.4275E-01, 
( 3.3105E+00, 
C G.0998E+01# 
( 7. 0620E+0 1 9 
C 1.3023E+03, 
C * * 
( * > 


7 

0) 

6.2084E+01) 
7.282GE+01) 
-2. 9069E+01 ) 
-1.2285E+02) 
-7. 1512E+01 ) 
-3.1662E+01) 
1.9G96E+01) 
* ) 

* ) 


( 1 . 0000E+G0* 
C 2. 0194E+0 1 9 
( 2.6518E-03, 
( 9. 4275E-01 > 
C 3.3105E+00* 
C 6.0998E+01* 
( 7.0620E+01, 
( 1 .3023E+03* 
C * » 
C * > 


8 

0) 

-6. 20S4E+01 ) 
-7 . 2826E+01 ) 
2. 9069E+0 1 ) 
1 . 2285E+02) 
7. 1512E+0 1 ) 
3. 1662E+01 ) 
-1 . 9696E+01 ) 
* ) 

* ) 
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Table 4.2 

Achievable Eigenvectors - Example 4.1 


ACH IEUABLE EIGEMUECTORS 
EIGENUALUES 

j: -1.4900 + J 2.3700 
2: -1.4900 - J 2.3700 
3: “• 0015 + J .0672 
4: -.0015 - J .0672 


EIGENUECTORS: (MAGNITUDE* PHASE (DEGREES) ) 


( 1 . QOOOE+OQ* 
( 2. 5129E+00 * 
( 1.9881E-02* 
( 8.g764E-01* 
( 9.0782E-03, 
( 5.2106E-02* 
( 2.5414E-02, 
( 1.4587E-01* 
C 1.5658E-01* 
( 4.4994E+00* 


( 1.0000E+00, 
( 3.2734E-01* 
( 2.39G7E+00, 
( 4 . 8700E+00 * 
( 2.3473E-02* 
( 5.1546E-02, 
( 1.5778E-03* 
( 3.4647E-03, 
( 1.356GE-01, 
( 3.7290E+00, 


1 

0) 

9.8272E+01) 
4. 5495E+01 ) 
-2. 3886E+01 ) 
2.7453E+01) 
-1 • 7034E+02) 
1 . 4961E+02) 
-4.8185E+01) 
-1 . 5405E+02) 
-1 . 5586E+02) 


3 

0 ) 

1.7157E+02) 
-1.7915E+02) 
8.0312E+01) 
2. 0170E+00 ) 
1 . 8022E+02) 
9. 3279E+01 ) 
-8.8518E+01) 
1 . 8008E+02) 
1 . 8057E+02) 


( 1.0000E+00, 
( 2. 5129E+00* 
( 1 . 9881E-02* 
( 8.97G4E-01* 
( 9. 0782E-03* 
( 5.210GE-02* 
C 2.5414E-02* 
( 1 . 4587E-01 * 
( 1 . 5658E-01 * 
( 4. 4994E+00* 


2 

0) 

-9. 8272E+01 ) 
-4.5495E+01) 
2. 388GE+01 ) 
-2* 7453E+01 ) 
1 . 7034E+02) 
-1.49G1E+02) 
4. 8185E+01 ) 
1 . 5405E+02) 
1 . 558GE+02 ) 


4 

0) 

-1.7157E+02) 
1 . 7S15E+02) 
-8.0312E+01) 
-2.0170E+00) 
-1.8022E+02) 
“9. 3279E+0 1 ) 
8 . 8518E+01 ) 
-1.8008E+02) 
-1 . 8057E+02 ) 


( 1 • 0000E+00* 

( 3.2734E-01, 

( 2.3967E+00* 

( 4. 8700E+00* 

( 2.3473E-02, 

( 5. 1546E-02* 

( 1.5778E-03* 

( 3.4G47E-03, 

( 1 . 3566E-01 * 

( 3.7290E+00, 


State vector (unsealed): x T = ( a , 0 r , u f , S r , ft , ) 
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Table 4.2, concluded 


EIGENUALUES 
5S -1,7577 + J 

BS -1,7577 - J 

7: -.5339 + J 

8: -.5339 - J 


8.6109 

8.6109 

21.3480 

21.3480 


EIGEhUECTORS: (MAGNITUDE* PHASE C DEGREES ) ) 


( 1.0000E+Q0, 
( 1.982BE+00, 
( 3.8565E-03, 
( 2.2560E-01, 
( 2.3537E+01, 
C 1 .OO35E-01* 
( 2.0685E+02, 
C 8.8194E-01, 
( 7. 1912E+00* 
( 2.2736E+02, 


5 

0) 

9.2299E+01) 
7.6557E+01) 
-9. 2379E+00 ) 
7. 0103E+01) 
3. 6410E+0 1 ) 
1 . 7164E+02) 
1.3795E+02) 
6. 1618E+0 1 ) 
-5.7208E+01) 


C 1.0000E+00, 
( 1 . 9826E+00 ► 
C 3.8565E-03, 
( 2.25B0E-01* 
( 2.3537E+01, 
C 1 . 0035E-01 » 
C 2.0685E+02, 
( 8.8194E-01, 
( 7 . 1912E+0Q* 
( 2.2736E+02, 


( 1.0000E+00* 
( 1.9760E+01, 
( 2.G355E-03, 
( 9.2533E-01, 
( 2.4966E+00, 
C 4.5728E+01* 
( 5.3314E+01* 
C 9.7650E+02, 
( 5.5269E+00, 
( 1.3441E+02. 


7 

0) 

6.7263E+01 ) 
7.5514E+01) 
-2.4169E+01) 
-1 . 3127E+02) 
-7.9773E+01) 
-3. 9837E+0 1 ) 
1 . 1659E+0 1 ) 
6. 5516E+00 ) 
-1 . 7060E+02) 


( 1.0000E+00. 
( 1.9760E+01* 
( 2.B355E-03* 
( 9.2533E-01, 
( 2 . 4966E+00# 
C 4. 5728E+0 1 p 
C 5. 3314E+0 1 * 
( 9 . 7650E+02* 
( 5 . 5269E+00 » 
C 1 . 3441E+02* 


6 

0) 

•9.2299E+01) 
7 . 6557E+01 ) 
9.2379E+00) 
*7. G103E+01 ) 
•3.B410E+01) 
■1.7164E+02) 
■1.3795E+02) 
■6.1S18E+01) 
5. 7208E+01 ) 


8 

0) 

G.7263E+01) 
•7 . 5514E+01 ) 
2. 4169E+01 ) 
1 . 3127E+02) 
7.9773E+01) 
3.9837E+01 ) 
■1.1B59E+01) 
•6.5516E+00) 
1 . 7060E+02 ) 
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(wf°) T = (1.0, 0.0) 

(w 2 °°) T = (0.0, 1.0) 


The control weighting matrix, R, can now be calculated by using Equation 
(4.11). The state weighting matrix, Q, can be calculated using Equations 
(4.13) through (4.16). These matrices are given in Table 4.3. 


Closed Loop System Analysis 

Four LQ control laws will be synthesized for the following values of con. 
trol weighting, p — ( 1.0 x 10 8 , 1.0 x 10 -2 , 1.0 x 10 -4 , 1.0 x 10 -8 ). 


Control synthesis for p — 1.0 x 10 8 . The LQ synthesis for this large value of p 
yields a closed-loop system whose eigenspace is close to that of the open-loop 
system. 

The system closed-loop eigenvalues are: 

X sp = - 1.348 ± j 2.193 = (w n = 2.574 (rad/s) , f= 0.524) 

X ph = - 0.0001 ± j 0.0526 = (w n = 0.053 (rad/s) , f= 0.002) 


\ = ~ 0.7264 ± j 8.758 = (w n = 8.788 (rad/s) , ?= 0.083) 
X ?2 = - 0.4564 ± j 21.351 = (w n = 21.36 (rad/s) , f= 0.021) 


X d = -9.0 
X cv = -10.0 


The eigenvectors of (A + BG) are given in Table 4.4. The gain matrix to 



Table 4.3 

State and Control Weighting Matrices 
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Table 4.4 

Closed-loop Eigenvectors - Example 4.1, p = 10 6 


LGEA RHQ=1 . 0E+06 


EIGENUALUES 
is -1.3479 + J 

2: -1.3479 - J 

3: -.0001 + J 

4: -.0001 - J 


2.1929 
2.1929 
.0526 
. 0526 


EIGENVECTORS: (MAGNITUDE* PHASE ( DEGREES ) ) 


( 1 . OOO0E+OO* 
( 2.2948E+00* 
( 2.1474E-02* 
( 8.9153E-01* 
( 2.8846E-01, 
( 5.1928E-02* 
C 7 . 4251E-01 * 
( 1 • 3366E-01 » 
( 8.0402E-08* 
( 9.6495E-08* 


( 1. 0000E+00, 
( 1.797GE-01* 
( 2. 1292E+00, 
C 3.4202E+00, 
( 2.5276E-01* 
( 5.2150E-02, 
( 1 . 3284E-02, 
( 2. 7409E-03, 
( 8. 6888E-05, 
( 2. 0050E-05* 


1 

0) 

9.5837E+01) 
4.5165E+01) 
-2.5741E+01) 
2.5619E+01) 
-1.7132E+02) 
1 • 4720E+02) 
-4.9747E+01) 
1.3879E+02) 
-1.4675E+02) 


3 

0) 

-1.9256E+02) 
-1 • 7933E+02) 
-2.8268E+02) 
3.7983E-01) 
-1.7984E+02) 
-2.6949E+02) 
-8.9716E+01) 
-2.7297E+02) 
-3.3516E+01) 


( 1 • OOOOE+OO* 
( 2.2948E+00* 
C 2.1474E-02* 
( 8.9153E-01* 
( 2.8846E-01* 
( 5. 1928E-02* 
( 7.4251E-01* 
( 1.3366E-01* 
( 8.0402E-08* 
( 9.6495E-08* 


C 1. OOOOE+OO* 
( 1.7976E-01* 
( 2. 12B2E+00* 
( 3. 4202E+00* 
( 2.5276E-01* 
( 5. 2150E-02* 
C 1.3284E-02* 
( 2.7409E-03* 
( 8.6888E-05, 
( 2.0050E-05* 


2 

0) 

-9.5837E+Q1) 
-4.5165E+01 ) 
2.5741E+01) 
-2.5619E+01) 
1.7132E+02) 
-1.4720E+02) 
4.9747E+01 ) 
-1.3879E+02) 
1.4675E+02) 


4 

0 ) 

1 . 9256E+02 ) 
1 • 7933E+02) 
2.8268E+02) 
-3.7983E-01) 
1 • 7984E+02) 
2.6949E+Q2) 
8. 9716E+01 ) 
2. 7297E+02) 
9.3516E+01 ) 


State vector (unsealed): X T = ( a , 0 T , u f , 6 r , & , & , ft , £2 , S e , S„ ) 
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Table 4.4, continued 


EIGENUALUES 

5: -.7264 + J 8.7584 
6: -.7264 - J 8.7584 
7: -.4564 + J 21.3506 
8: -.4564 - J 21.3506 


E IGENUECT0R5 : ( MAGN ITUDE* PHASE ( DEGREES ) ) 


( 1.0000E+00, 

( 7.7404E+00, 

( G.2077E-03, 

C 8.8074E-01, 
( 8.8783E+00, 

( 2.43S0E-O1# 

C 7.8027E+01, 
( 2.1435E+00, 

( 2. 7734E-06* 

( 8. 8981E-07 9 


5 

0 ) 

-2.9187E+02) 
-2.8869E+02) 
-2.6606E+01) 
-2.3329E+02) 
-2. 191GE+02) 
-1.3855E+02) 
-1.2442E+02) 
-3.3344E+02) 
-3. 4322E+02) 


( 1.0000E+00, 
C 7.7404E+00, 
( G.2077E-03, 
C 8.8074E-O1* 
C 8.8783E+00, 
( 2.4390E-01, 
C 7.8027E+01* 
( 2.1435E+00, 
( 2.7734E-06* 
( 8.8981E-07 9 


( 1.0000E+00, 

( 2.0127E+01, 

C 2.6356E-03, 
C 9.4247E-01, 
C 3.3014E+00, 
C 6.0949E+01, 
( 7. 0502E+01 » 

( 1 . 3016E+03* 

C 2.4073E-07, 
C 3.1332E-07, 


7 

0 ) 

G.2180E+01) 
7.2932E+01) 
-2.9045E+01) 
2.3713E+02) 
-7 . 151GE+01 ) 
-3.1G47E+01) 
1 • 9709E+01 ) 
-3.1940E-MD1) 
2. 0149E+02) 


( 1 • 0000E+00f 
( 2. 0127E+01* 
( 2.6356E-03, 
( 9.4247E-01, 
( 3.3014E+00, 
( 6. 0949E+01 t 
( 7 . 0502E+01 * 
( 1.3016E+03, 
( 2.4073E-07, 
( 3.1332E-0 7 , 


6 

0 ) 

2. 9187E+02) 
2.8869E+02) 
2. 660GE+01 ) 
2.3329E+02) 
2. 1916E+02) 
1.3855E+02) 
1 . 2442E+02) 
3. 3344E+02) 
3. 4322E+02) 


8 

0) 

6.2180E+01 ) 
•7.2932E+01) 
2. 9045E+01 ) 
■2.3713E+02) 
7 . 1516E+01) 
3. 1G47E+01 ) 
1.9709E+01) 
3. 1940E+01 ) 
■2. 0149E+02) 
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Table 4.4, concluded 


EIGENVALUES 

9t -9. 0000 + J 0 

10: -10.0000 + J 0 


EIGENVECTORS: (MAGNITUDE* PHASE ( DEGREES ) ) 


10 

0) 

1 . 8000E+02) 
0) 
0) 
0) 
0) 

1.80Q0E+02) 

1.S000E+02) 

0) 

0) 


( 1 . QQ0QE+00* 
( 9.2712E+00, 
( 6.8055E-03, 
( 1 . 0301E+00* 
( 2. 27S1E+D0* 
( 1 . 5883E-01 * 
( 2.Q4&5E+01, 
( L.4295E+G0, 
( 5. 4943E+00* 
( 5.2G24E-08, 


9 

0) 

-1.8000E+02) 

0) 

0) 

-1.8000E+02) 

0) 

0) 

-1.8000E+02) 

-1.8000E+02) 

0) 


( 1 . 0000E+00* 
( 3.142GE+00, 
( 5.7315E-03, 
( 3.142SE-01, 
( 2.7597E+00, 
( 1.1082E-01. 
( 2. 7597E+01 * 
( 1 . 1082E+00* 
( 7.302BE-06* 
( 1.9544E+02, 
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obtain this eigenspace in the closed-loop system and the closed-loop system 
matrices for this synthesis are given in Appendix A.4. The closed-loop system 
residue magnitudes, for u p = impulse, for the short period, phugoid, first elas- 
tic, and second elastic modes for each output are shown in Figure 4.1. Again 
note that the impulse residues are normalized for plotting such that the sum of 
the short period, phugoid, first elastic, and second elastic modes residue magni- 
tudes in each response is unity. The output time responses due to a unit step 
in Up, with zero initial conditions, are shown in Figure 4.2. 


Control synthesis for p = 1.0 x 10 2 . The closed-loop system eigenvalues for 
this value of control weighting are: 

X 12 = “ 2.264 ± j 1.137 = (w n = 2.534 (rad/s) , 0.894) 

X 3>4 = - 0.0490 ± j 0.1464 = (w n = 0.154 (rad/s) , f= 0.317) 

X 66 = - 7.385 ± j 10.87 = (w n = 13.14 (rad/s) , f= 0.562) 

X 78 = - 0.5396 ± j 21.45 = (w n = 21.46 (rad/s) , f= 0.025) 


X 9 = - 10.05 
X 10 = - 13.22 

The eigenvectors of (A + BG) are given in Table 4.5. The gain matrix and 
the closed-loop system matrices for this value of p are given in Appendix A.4. 
The normalized closed-loop system impulse residue magnitudes of the short 
period, phugoid, first elastic, and second elastic modes for each output are 
shown in Figure 4.3. The output time responses due to a unit step in n p , with 
zero initial conditions, are shown in F igure 4.4. 



MRGNITUDE MR0NITUDE MRGNITUDE MAGNITUDE 
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IMPULSE RESIDUE MRGNITUDES 
(DUE T0 PILOT INPUTS) 

LQER RH0=1.QE+Q6 

GfiMMfl NZP 



MODES MODES 

THETfl-R THETfl-R DOT 



MOOES 


MODES 


U-VELOCITY 



E2 El SP PH 
MODES 


Figure 4.1 

Impulse Residue Magnitudes - Example 4.1, p = 10 8 
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Figure 4.2 

Output Time Responses - Example 4.1, p = 10 6 


10 2 . 

TIME (S) 

Figure 4.2, continued 
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Figure 4.2, concluded 
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Table 4.5 

Closed-loop Eigenvectors - Example 4.1, p = 10“ 2 


LQEA RH0=1 • 0E~02 


EIGENUALUES 
is “2.2642 + J 

2: -2.2642 - J 

3: -.0490 + J 

4: -.0490 - J 


1 . 1367 
1.1367 
.1464 
.1464 


EIGEMUECTORS: (MAGNITUDE* PHASE (DEGREES) ) 


( 1 . 0000E+00* 
( 1.6980E+00, 
( 1.3235E-02* 
( 6.7022E-01, 
( 2.6766E-01* 
( 3.4672E-02* 
( 6.7812E-01, 
( 8.7844E-02* 
( 5. 4025E“0i* 
C 1.5166E+00* 


( 1.0000E+00, 
C 2. 4521E+00* 
( 3. 4442E+00* 
C 1 . 5885E+01 * 
( 7.9536E“01> 
( 3.7454E“02* 
C 1.2277E-01* 
( 5. 7815E-03* 
( 1 • 4265E+00* 
( 1 . 8642E+00* 


1 

0) 

1 . 3474E+02) 
1.8329E+01) 
“1 . 8603E+01 ) 
1 . 5226E+02) 
1.9613E+02) 
“5.4400E+01) 
-1.0530E+01) 
1.6477E+02) 
-1.4788E+02) 


3 

0) 

“5. 1894E+01 ) 
-7 . 8686E+01 ) 
-1.6038E+02) 
1.2833E+02) 
-1.3425E+02) 
-1.2318E+02) 
“2. 5756E+01 ) 
1 . 3614E+02) 
“1 • 2374E+02) 


( 1 . 0000E+00* 
( 1 • 6980E+00* 
C 1 . 9235E-02* 
( 6. 7022E-01* 
( 2.6766E-01* 
( 3.46722-02* 
( 6.7812E“01 * 
( 8.7844E-02* 
( 5.4025E-01, 
( 1.5166E+00* 


( 1.000QE+00. 
( 2.4521E+00, 
( 3.4442E+00* 
( 1.5885E+01, 
( 7.9536E-01* 
( 3.7454E-02* 
( 1.2277E-01* 
( 5.7815E-03, 
( 1 . 4265E+00* 
( 1.8B42E+00, 


2 

0) 

-1.3474E+02) 
-1.8329E+G1) 
1.8603E+01) 
-1.5226E+02) 
-1.9B13E+02) 
5. 4400E+01 ) 
1 . 0530E+01 ) 
-1.6477E+02) 
1.4788E+02) 


4 

0 ) 

5. 1894E+01 ) 
7.8686E+01) 
1 . 6038E+02) 
“1 . 2833E+02) 
1 • 3425E+02) 
1.2318E+02) 
2. 5756E+01 ) 
“1 .3614E+02) 
1.2374E+02) 


State vector (unsealed): x T = ( a , 0 T , u f , 0 r , ^ > kx » £2 » Se , ) 
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Table 4.5, continued 


EIGENUALUES 
5: -7.3843 

G: -7.3843 

7: -.5336 

8s 533G 


+ J 10.863G 

- J 10.8636 

+ J 21.4516 

- J 21.4516 


EIGENUECTORS: (MAGNITUDE* PHASE ( DEGREES ) ) 


( 1.0000E+00* 
C 1.3601E+01* 
( 4 .5932E-03* 
( 1.0350E+00, 
( 4.4418E+00* 
( 5. 7334E-01* 
( 5.837GE+01, 
( 7.5343E+00, 
( 1 • 0053E+01 > 
( 1 . 2467E+01 * 


C l.OOOOE+OO 
( 2.0132E+Q1 
( 2.6205E-03 
( 3.3821E-01 
( 2. 4657E+00 
( 3. 9175E+01 
( 5.2909E+01 
( 8 . 40S2E+02 
( 3.3608E+00 
( 1.3773E+0I 


5 

0 ) 

1 • 0380E+02) 
4.4263E+01) 
3 . 39S1E+02) 
1 . 3739E+02) 
2. 4047E+02) 
2 . 6158E+02) 
4 . 6582E+00 ) 
4. 9394E+01 ) 
1.291GE+01) 


7 

0) 

6 . 3106E+01 ) 
7. 3138E+01 ) 
-2. 8335E+0 1 ) 
-1.5S94E+02) 
-7.0820E+01) 
-6. 5499E+01 ) 
2.0621E+01) 
-3. 540QE+01 ) 
-1.6233E+02) 


( 1.0000E+00* 
( 1 . 3601E+01 * 
( 4.5932E-03* 
( 1 . 0350E+00* 
( 4. 4418E+00* 
( 5.7334E-01* 
C 5.8370E+01* 
( 7. 5343E+00 * 
( 1 . 0053E+01 * 
( 1 . 24G7E+01 * 


( 1.0000E+00* 
( 2. 0132E+01 * 
C 2.6205E-03* 
( 9 . 3821E-Q1 * 
( 2. 4657E+00* 
( 3. 9175E+01* 
( 5. 2909E+01 * 
( 8 . 4062E+02* 
( 9. 3608E+00* 
( 1.3773E+01* 


6 

0) 

1.0380E+02) 
■4.4263E+01) 
•3.39S1E+02) 
•i .3739E+02) 
2.4047E+0R) 
2.6158E+02) 
4.G582E+00 ) 
4. 9394E+01 ) 
1 . 2916E+0i ) 


8 

0) 

G.3106E+01) 
7 . 3138E+01 ) 
2. 8335E+01 ) 
1 . 5694E+02) 
7 . 0820E+01 ) 
6. 5499E+0 1 ) 
•2.0621E+01) 
3.5400E+01 ) 
1.6233E+02) 
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Table 4.5, concluded 


EIGENUALUES 

9: -10.0486 + J 0 

10: -13.2191 + J 0 


EIGENVECTORS: (MAGNITUDE, PHASE C DEGREES ) ) 


C 1 . OOOOE+OO > 
( 1.0780E+01, 
( G. 2370E-03, 
( 1.0728E+00, 
C 2.6042E+00, 
( 2. 1026E“0 1 r 
C 2.61G8E+G1, 
( 2. 1 128E+00, 
( 7.0229E+00, 
( 7. 9220E-01 , 


9 

0) 

1.8000E+025 

0) 

0) 

1.8000E+02) 

0) 

0) 

1.3000E+02) 
i . 8000E+02) 
1.8000E+Q2) 


C 1. 0000E+00, 
C 1.6138E+01, 
( 5.I1G6E-Q3, 
( 1 . 2208E+00, 
C 3. 9217E+00 , 
( 4.0887E-01, 
( 5. 1841E+01 * 
C 5. 4049E+Q0, 
C 1.4050E+01, 
C 2.G095E+01, 


10 

0) 

1 • 8000E+02) 
0) 
0 ) 

1.8000E+02) 

0) 

0) 

1 , 8000E+02) 
1 . 8000E+02) 
1.8000E+02) 



MAGNITUDE MAGNITUDE MAGNITUDE MAGNITUDE 
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IMPULSE RESIDUE MAGNITUDES 
(DUE TO PILOT INPUTS) 
LQEA RH0=1 .0E-02 
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Figure 4.3 

Impulse Residue Magnitudes - Example 4.1 , p = 10 -2 
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Figure 4.4, continued 



THETR-T DOT ( RRD/S ) (XIO -*) THETR-T (RRD) (XIO “*) 



Figure 4.4, continued 


NZP (G'S) 
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Figure 1.4, concluded 
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Control synthesis for p - 1.0 x 10 4 . The system closed-loop eigenvalues, the 
eigenvalues of (A + BG) are: 

Xj 2 = - 1.696 ± j 2.064 = (w n = 2.672 (rad/s) , 0.635) 

X 3 4 = - 0.0445 ± j 0.0971 = (w n = 0.107 (rad/s) , 0.417) 

X 5 6 = — 4.079 ± j 8.355 = (u; n = 9.297 (rad/s) , 0.439) 

X 7j8 = - 0.5959 ± j 21.46 = (w n = 21.47 (rad/s) , $•= 0.028) 


Xg = - 54.11 
Xjq — 98.55 

The eigenvectors of (A + BG) are given in Table 4.6. The gain matrix and 
the closed-loop system matrices for this value of p are given in Appendix A.4. 
The normalized closed-loop system impulse residue magnitudes of the short 
period, phugoid, first elastic, and second elastic modes for each output are 
shown in Figure 4.5. The output time responses due to a unit step in u p , with 
zero initial conditions, are shown in Figure 4.6. 


Control synthesis for p — 1.0 x 10 6 . The system closed-loop eigenvalues, the 
eigenvalues of (A + BG) are: 

X sp = - 1.493 ± j 2.365 = (w n = 2.797 (rad/s) , £= 0.534) 

X ph = - 0.0075 ± j 0.0679 = (w n = 0.068 (rad/s) , 0.110) 

X^= - 1.807 ± j 8.608 = (w n = 8.796 (rad/s) , 0.205) 
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Table 4.6 

Closed-loop Eigenvectors - Example 4.1, p = 10 4 


LQEA RHQ=1 . 0E-04 


EIGENUALUE5 
l: -1.6957 + J 

E2 -1.6957 - J 

32 -.0445 + J 

4: -.0445 - J 


2.0645 

2.0645 

.0971 

.0971 


EIGEMUECTOR52 (MAGNITUDE, PHASE (DEGREES) ) 


( 1 . 0000E+00 » 

( 2.2485E+00, 

( 2.0133E-02, 

( 8.41G4E-01, 

( 9. 2655E-02* 

( 4 . 7214E-02, 

( 2.4754E-01, 

( 1.2614E-01, 

( 2. 5761E-0 1 * 

( 3.7742E+00, 


( 1 . 0000E+00, 
( 9. 0230E-01 , 
( 2.6739E+00, 
( 8.4477E+00, 
( 2. 7688E-01 * 
( 5.09G7E-02, 
( 2.9574E-02, 
( 5.4438E-03, 
( 4 . 9927E-01 » 
( 3.6476E+00, 


1 

0 ) 

1 . 0495E+02) 
3. 8333E+01 ) 
-2.4447E+01) 
1.4364E+02) 
-1.6760E+02) 
-8. G9G0E+01 ) 
-3. 8206E+01 ) 
-1.8363E+Q2) 
-1 . 5251E+02) 


3 

0) 

2.5597E+02) 
2. 1640E+02) 
1 . 4134E+02) 
9.3930E+01) 
1.9449E+02) 
2. 0857E+02) 
3. 0913E+02) 
1 . 1451E+02) 
1 . 9670E+02) 


( 1.0000E+00, 
( 2.2485E+00, 
( 2. 0139E-02, 
( 8.41G4E-01, 
( 9.2655E-02, 
( 4.7214E-02, 
C 2.4754E-01, 
( 1 . 2614E-0 1 t 
( 2.5761E-01, 
C 3.7742E+00, 


( 1 . 0000E+00 , 
( 9.0230E-01, 
( 2. G739E+00, 
( 8.4477E+00, 
( 2.7G88E-01, 
( 5.0967E-02, 
( 2. 9574E-02, 
( 5.4438E-03, 
( 4.9927E-01, 
( 3. 6476E+00, 


2 

0) 

-1.0495E+02) 
-3.8333E+01) 
2.4447E+01) 
-I . 4364E+02) 
1.G760E+02) 
8. G9GOE+01 ) 
3.8206E+01 ) 
1.8363E+02) 
1 . 5251E+02) 


4 

0) 

-2. 5597E+02) 
-2.1640E+02) 
“1 . 4134E+02) 
-9. 3930E+01 ) 
-1 . 9449E+02) 
-2.0857E+02) 
-3. 0913E+02) 
-1.1451E+02) 
-1 . 9G70E+02) 


State vector (unsealed): X T = ( a , S T , u f , 0 T , & , &> > ^2 , > <*cv ) 
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Table 4.6, continued 


EIGENUALUES 

5: “4*0793 + J 8.3546 
6: -4.0793 - J 8.3546 
7 s -.5959 + J 21.4568 
8: -.5959 - J 21. 4568 


EIGEMUECTORS: (MAGNITUDE, PHASE (DEGREES) ) 


( 1 • 0000E+00, 
( 7.9835E+00, 
( 5.8743E-03, 
( 8.58S9E-01, 
( 6. 0735E+00, 
( 2.4303E-01, 
( 5. G467E+01 , 
( 2. 2595E+00* 
( 4. 9248E+00, 
( 4.86G0E+01, 


5 

0) 

9. 7143E+01 ) 
5.4223E+01 ) 
-1 • 8882E+01 ) 
1 . 0012E+02) 
-1 • 5274E+02) 
-1.438GE+02) 
-3. G716E+01 ) 
4.3704E+01) 
-2.3108E+01) 


( 1 . 0000E+00 
( 7.9835E+00 
( 5.8743E-03 
( 8.5869E-01 
( 6. 0735E+00 
( 2.4303E-01 
( 5. G467E+01 
( 2. 2595E+00 
( 4. 9248E+00 
( 4. 86G0E+01: 


( 1.0000E+00, 
( 1.998GE+01, 
( 2.6208E-03, 
( 9.3108E-01, 
( 2. 2055E+00, 
( 3. 3925E+01 , 
( 4. 7342E+01 , 
( 7 • 2821E+02, 
( 1 . 0542E+01 , 
( 7. 5550E+01 , 


7 

0) 

G.544GE+01) 
7 . 4240E+01 ) 
-2. G144E+01 ) 
1.9844E+02) 
-7.4597E+01) 
-G.9968E+01) 
1 . G394E+01 ) 
-2.7466E+01) 
1.9043E+02) 


( 1 • 0000E+00, 
( 1.9986E+01, 
( 2.6208E-03, 
( 9.3108E-01, 
( 2. 2055E+00, 
( 3.3925E+01, 
( 4. 7342E+01 , 
( 7 . 2821E+02, 
( 1.0542E+01, 
( 7 . 5550E+01, 


G 

0) 

-9.7143E+01) 
■5.4223E+01) 
1.8882E+01) 
■1.0012E+02) 
1.5274E+02) 
1.438GE+02) 
3.G71GE+01) 
4.3704E+01 ) 
2.3108E+01 ) 


8 

0) 

G. 5446E+01 ) 
■7 . 4240E+01 ) 
2.6144E+01) 
1.9844E+02) 
7.4597E+01 ) 
6. 9968E+01 ) 
•1.G994E+01) 
2.74G6E+01 ) 
1.9043E+02) 



106 


Table 4.6, concluded 


EIGENUALUES 

3: -54.1139 + J 0 

10: -93.550? + J 0 


EIGEMUECTORS: (MAGNITUDE* PHASE (DEGREES) ) 


10 

0) 

0) 

0) 

1.8000E+02) 

0) 

1 . 8000E+02) 

1 . 800GE+02) 
0) 
0) 
0) 


C l.OOOOE+OOt 
C 6.234GE+02* 
( 6.6947E-03, 

( i . 1520E+01 * 
( 4.588GE+01* 

( j. . 3822E+01 * 
( 2.4831E+03, 

( 7. 4798E+02* 

( 2. 1747E+03* 

( 1.8788E+03. 


9 

0) 

0) 

-1.8000E+02) 

-1.8000E+02) 

0) 

-1.8000E+02) 

-1.8000E+02) 

0 ) 

0 ) 

0) 


( 1.0000E+00, 
( 5.731 1E+01* 
( G.8349E-05, 
( 5.8154E-01* 
( 3.8181E+00* 
( 8.8285E-01, 
C 3.7628E+02, 
( 8 . 7006E+G 1 9 
( 4.5G49E+02* 
( 4.0420E+03* 



MAGNITUDE MAGNITUDE MAGNITUDE MAGNITUDE 
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IMPULSE RESIDUE MRGNITUDES 
(DUE TO PILOT INPUTS) 
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Figure 4.5 

Impulse Residue Magnitudes - Example 4.1, p — 10 -4 
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Figure 4.6, continued 


THETR-T DOT ( RRD/S ) CXIO - 1 ) THETR-T (RflD) (X10 - 1 ) 



Figure 4.6, continued 




NZP (G'S) 
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Figure 4.6, concluded 
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X^= - 0.5351 ± j 21.35 = (w n = 21.36 (rad/s) , f= 0.025) 


X el = - 500.5 
X cv = - 999.9 

The eigenvectors of (A + BG) are given in Table 4.7. The gain matrix and 
the closed-loop system matrices for this value of p are given in Appendix A.4. 
The normalized closed-loop system impulse residue magnitudes of the short 
period, phugoid, first elastic, and second elastic modes for each output are 
shown in Figure 4.7. The output time responses due to a unit step in u p , with 
zero initial conditions, are shown in Figure 4.8. 


Comments on Example 4.1 

This example illustrates the use of the LQEA control synthesis technique. 
LQEA yields a family of control laws as a function of LQ cost function control 
weighting. Each of these control laws will be discussed separately. 

The synthesis for p = 1.0 x 10 6 , as is expected, yields a closed-loop system 
whose eigenspace is close to that of the open-loop system. This can be seen by 
comparing the impulse residue magnitudes and time responses obtained from 
this synthesis with those of the open-loop system presented in Chapter III. 

The control synthesis for p = 1.0 x 1(T 2 and for p - 1.0 x 10 -4 yields 
interesting results. By comparing the closed-loop system eigenvectors for each 
value of control weighting, one can see how the closed-loop system eigenvectors 
change as they tend towards their asymptotic values. For these values of con- 
trol weighting, the closed-loop eigenspace is not close to the desired eigenspace. 
The impulse residue magnitudes and time responses show that the first elastic 
mode is contributing significantly to the rigid-body pitch and rigid-body pitch 
rate outputs. The control laws for these two values of control weighting are 
actually increasing the contribution of the first elastic mode in these two out- 
puts. This is in opposition to the desired synthesis result. 

A control weighting of p = 1.0 x 10~ 6 is required to achieve the desired 
synthesis objectives. The short period, first elastic, and second elastic eigen- 
values are close to the specified desired finite eigenvalues and the eigenvectors 
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Table 4.7 

Closed-loop Eigenvectors - Example 4.1, p 


LQEA RH0=1 • 0E-06 


EIGENVALUES 
is -1.4930 + J 

SS -1.4930 - J 

3: -.0075 + J 

4s -.0075 - J 


2.3655 

2.3655 

.0679 

.0679 


EIGENVECTORS: (MAGNITUDE* PHASE (DEGREES) ) 


1 

( 1 • OOOOE+OO* 0) 
( 2.5088E+00, -2.6164E+02) 
( 1.9887E-02, -3. 1461E+02) 
( 8. 9688E-01 * -2.3899E+01 ) 
( 8.4733E-03* -3.2331E+02) 
( 5. 2032E-02* -1 . 7030E+02) 
( 2.3702E-02, -2.0105E+02) 
C 1.4555E-01* -4.8041E+01) 
( 1.5756E-01* -1 . 5479E+02) 
( 4.4881E+00* -1 • 5581E+02) 


( 1 . OOOOE+OO* 
( 3. 3618E-01 » 
( 2.4097E+00, 
( 4.9237E+00* 
( 3.4129E-02* 
( 5.1619E-02* 
( 2.3303E-03* 
( 3.5245E-03* 
( 1.3758E-01, 
C 3.7319E+00* 


3 

0) 

1.8462E+02) 
1 • 8378E+02) 
8. 8301E+0I ) 
4.2316E+01) 
1 . 8134E+02) 
1 . 3864E+02) 
2. 7767E+02) 
i .6462E+02) 
1 . 8184E+02) 


( 1 • 0000E+00* 
( 2.5088E+00* 
( 1.9887E-02* 
( 8.9688E-01, 
( 8.4733E-03, 
( 5.2032E-02* 
( 2.3702E-02, 
( 1.4555E-01, 
C 1.5756E-01, 
( 4. 4881E+00* 


( 1 • 0000E+00* 
( 3.3618E-01* 
( 2. 4097E+00* 
( 4.9237E+00* 
( 3.4129E-02* 
( 5.1619E-02* 
( 2.3303E-03* 
( 3.5245E-03, 
( 1.3758E-01, 
( 3.7319E+00, 


State vector (unsealed): x T = ( a , 6 T , u f , 6 r , ft , > 


2 

0 ) 

2.6164E+02) 
3. 1461E+02) 
2.3899E+01 ) 
3.2331E+02) 
1.7030E+02) 
2. Q105E+02) 
4. 8041E+01 ) 
1.5479E+02) 
1.5581E+02) 


4 

0) 

■1.8462E+02) 
1 • 8378E+02) 
8. 8301E+01 ) 
4. 2316E+01 ) 
1 . 8134E+02) 
1 . 3864E+02) 
2.7767E+02) 
1 • 6462E+02) 
1 . 8184E+02) 


£2 i i ^cr ) 
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Table 4.7, continued 


EIGENUALUE5 


5: 

-1.80B5 + J 

8.6082 

B: 

-1.80B5 - J 

8.6082 

7: 

-.5351 + J 

21.3497 

8: 

-.5351 - J 

21.3497 


EIGENVECTORS: (MAGNITUDE* PHASE (DEGREES) ) 


( 1 • 0000E+00* 

( 2.4172E+00* 

( 4.5304E-03* 

( 2.7482E-01* 

( 2.2073E+01* 

( 8it)5G8E“02* 

( 1.3415E+02* 

( 7. 08B5E-01 * 

( 6.8016E+00, 

( 2. 1259E+02* 


5 

0 ) 

-2.713IE+02) 
“2.8511E+02) 
-1.3167E+01) 
-2.8878E+02) 
-3.1145E+02) 
-1 .8692E+02) 
-2. 0960E+02) 
-2.9977E+02) 
-5.6019E+01 ) 


( 1.0000E+00, 
( 2.4172E+00* 
( 4.0304E-03* 
( 2.7482E-01* 
( 2.2073E+01 > 
( 8.0568E-02* 
( 1.9415E+02* 
( 7 . 0865E-01 » 
( 6.8016E+00* 
( 2. 1259E+02* 


( l.OOOOE+OOt 
( 1 .9764E+01 * 

( 2. B353E-03* 

( 9.2545E-01 > 

( 2.483'BE+OO, 

C 4.5460E+01* 
( 5.3044E+01* 

( 9.7087E+02* 

( 5.S99BE+00* 

( 1.3312E+02, 


7 

0) 

G.7223E+01) 
7.5487E+01 ) 
“2. 4212E+01 ) 
-I • 3184E+02) 
-7.9G91E+01) 
-4. 0409E+01 ) 
1.1745E+01) 
5.2243E+00) 
-1.7059E+02) 


( 1 . OOOOE+OO* 
( 1.9764E+01* 
( 2.B353E-03* 
( 9.2545E-01* 
( 2.4838E+00* 
( 4. 54GOE+01 » 
( 5. 3044E+01 * 
( 9. 7087E+02* 
( 5.5996E+00, 
( 1.3312E+02. 


B 

0 ) 

2.7131E+02) 
2.8511E+02) 
1 . 3167E+01 ) 
2.8878E+02) 
3. 1145E+02) 
1 • 8692E+02) 
2.09G0E+02) 
2. 9977E+02) 
5.G019E+01) 


8 

0) 

•6.7223E+01 ) 
7.5487E+01 ) 
2.4212E+01 ) 
1.3184E+02) 
7.9691E+01 ) 
4. 0409E+01 ) 
•1.1745E+01) 
’5.2243E+00 ) 
1.7059E+02) 
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Table 4.7, concluded 


EIGENUALUES 

9: -500.4525 + J 0 

10S -939.8590 + J 0 


E I GENUECTORS : C M AGN I TUDE , PHASE (DEGREES ) ) 


( 1.0000E+00f 
( 5.7613E+01* 
( 4.4385E-05, 
( 1 . 1512E-0 1 f 
( 4.52G6E-01, 
( 1.6190E-01, 
( 2.2G54E+02, 
( 8. 1025E+01* 
C 1.9097E+03, 
C 1.9G54E+02, 


9 

0) 

0) 

0) 

1.8000E+02) 

0) 

1.8000E+02) 

1.8000E+02) 

0) 

0) 

0) 


( 1 . 0000E+Q0? 
( 2. 5169E+01 » 
( 2. 7197E-05* 
( 2.5172E-02, 
( 3. 9915E-01 » 
( 2.8332E-03, 
( 3.9909E+02, 
( 2.8328E+00, 
( 1.3820E+03* 
( 1 . 3023E+05* 


10 

0) 

1 , 8000E+02) 
0) 
0) 
0) 

1.8000E+Q2) 

1 . 8000E+02) 
0) 
0) 
0) 
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IMPULSE RESIDUE MRGNITUDES 
(DUE TO PILOT INPUTS) 
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Figure 4.7 

Impulse Residue Magnitudes - Example 4.1, p — 10 -6 
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Figure 4.8 

Output Time Responses - Example 4.1, p = 1(T 6 
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Figure 4.8, concluded 
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associated with these modes are close to the achievable eigenspace. This is 
reflected in the closed-loop system impulse residues and time responses. This 
control law is decreasing the contribution of the elastic modes to the rigid-body 
dynamics. Also, in this case, one can note that the augmented actuator eigen- 
values have achieved the specified two-to-one eigenvalue bandwidth ratio, and 
their associated eigenvectors have tended towards uncoupling the actuator 
dynamics. 

In this case, very large actuator bandwidths are required. The resulting 
control vane actuator bandwidth is 1000 (rad/s) and the elevator actuator 
bandwidth is 500 (rad/s). Though not possible, if the required actuator 
bandwidth could be achieved, it is clear that filtering would be required to 
prevent exciting higher order, unmodeled elastic modes. 

It should be noted that a low value of control weighting was requirec 
obtain satisfactory results. By comparing the 0 r time responses for each vaiv 
of control weighting, one can see that the amplitude of these responses tends i 
decrease as p decreases. It appears that as p is decreased, the pilot’s input 
having less effect upon the system’s dynamic response. This result will 
explored in a later section. 

A different choice of the desired eigenspace can lead to significantly 
different results. This is illustrated in the next example. 


Example 4.2 

This example illustrates how a change in one mode of the desired eigen- 
space can significantly affect the synthesis results. For this example, the 
desired infinite eigenvalues and eigenvectors are chosen to be the same values 
as in Example 4.1. The desired finite eigenvectors and eigenvalues are also 
chosen to be the same as in Example 4.1, with the exception of the second elas- 
tic mode. 

The second elastic eigenvalue is chosen to be placed at a location such 
that the mode’s natural frequency is the same as the open-loop natural fre- 
quency and the mode has increased damping. This value was chosen to be: 

X& = -4.271 ± j 20.92 = (w n = 21.35 (rad/s) , f = 0.200) 

The desired eigenvector for the second elastic eigenvector is chosen using the 
same strategy as used in Example 3.1 and 3.2. The elements of the eigenvector 
that are associated with the (scaled) rigid-body states are chosen to have small 
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magnitudes. The element associated with the state £ 2 * s chosen to be unit 
magnitude as a reference value. The chosen desired second elastic mode’s 
eigenvector is given in Table 4.8. The achievable second elastic eigenvalue is 
given in Table 4.9. The achievable eigenvectors for the short period, phugoid, 
and first elastic modes are identical to those given in Table 4.2. 

Cpptro| laws were synthesized for various values of control weighting, p. 
At a control weighting of p = 1.0 x 10 -6 , the LQEA synthesis yields a closed- 
loop system eigenspace that is close to the desired eigenspace. The eigenvalues 
for this SFnffrpsis arp: 

X sp = -1.495 ± j 2.366 =(a; n = 2.799 (rad/s) , = 0.534) 


X ph = -0.0075 ± j 0.0680 =(w n = 0.068 (rad/s) , £ = 0.109) 


X Cl = -1.774 ± j 8.638 =(w n = 8.818 (rad/s) , f = 0.201) 
X fz = -6.401 ± j 19.02 =(w n = 20.07 (rad/s) , f = 0.319) 

X el = -595.9 


X cv = -945.9 

The normalized closed-loop system impulse residue magnitudes for each output 
are shown in Figure 4.9. 


Comments on Example 4>2 

In this example, the eigenvector selection strategy that worked well with 
the short period, phugoid, and first elastic modes has yielded an achievable 
second elastic eigenvector, calculated from */| 2 and X| 2 , that is a “mixed mode”. 
As can be seen in Table 4.9, in this “second elastic” eigenvector the element 
associated with the element has a large magnitude. 

This control law yields closed-loop system eigenvalues that are close to the 
desired locations, but unacceptable impulse residue magnitudes. By comparing 
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Table 4.8 

Desired Second Elastic Eigenvector - Example 4.2 


DESIRED EIGENVECTORS 
EIGENVALUES 

i: -4.2709 + J 20.9230 

2: -4.2709 - J 20.9230 


EIGENVECTORS: (MAGNITUDE* PHASE ( DEGREES ) ) 




1 


2 

( 

1.0000E-05* 

0) 

( 1.0000E-05, 

0) 

( 

1.0000E-04* 

6. 0000E+01 ) 

( 1.0000E-04* 

-6. 0000E+01 ) 

( 

1.0000E-08* 

7. 0000E+01 ) 

( 1.0000E-08, 

-7 . OOOOE+Oi ) 

( 

1 . 0000E-05? 

-3. 0000E+01 ) 

( 1.0000E-05, 

3. 0000E+01 ) 

( 

* 9 

* ) 

( * » 

* ) 

( 

* 9 

* ) 

( * * 

* ) 

( 

* 9 

* ) 

( * * 

* ) 

( 

1. 0000E+00* 

2. 0000E+01 ) 

( 1. OOOOE+OO* 

-2. 0000E+01 ) 

( 

* , 

* ) 

( * 9 

* ) 

( 

* , 

* ) 

( * * 

* ) 


Table 4.9 

Achievable Second Elastic Eigenvector - Example 4.2 


ACHIEVABLE EIGENVECTORS 
EIGENVALUES 

l: -4.2709 + J 20.9230 

2: -4.2709 - J 20.9230 


EIGENVECTORS: (MAGNITUDE* PHASE (DEGREES ) ) 


( 1 . 0000E+00* 
( 8.5440E-02* 
( 1.2283E-03* 
( 4. 0010E-03* 
( 1 . 5204E+01 * 
( 2.7229E+00* 
( 3. 24G7E+02* 
( 5.8146E+01* 
( 3 . 0235E+01 * 
( 1.G273E+03* 


1 

0 ) 

-5. 2034E+01 ) 
-2.8160E+02) 
-1.5357E+02) 
-2.8515E+02) 
-1.9847E+02) 
-1.8361E+02) 
-9.6934E+01) 
-8. 387GE+01 ) 
-7.9623E+01) 


( 1 . OOOOE+OO* 
( 8.5440E-02* 
( 1.2283E-03* 
( 4.0010E-03* 
( 1.5204E+01, 
( 2.7229E+00, 
( 3.24G7E+02* 
( 5.8146E+01, 
( 3.0235E+01* 
( 1.6273E+03* 


2 

0 ) 

5.2034E+01) 
2.81E0E+02) 
1.5357E+02) 
2.8515E+02) 
i . 9847E+02) 
1.8361E+02) 
9. G934E+01 ) 
8.387GE+01) 
7.9G23E+01) 
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Figure 4.9 

Impulse Residue Magnitudes - Example 4.2 
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the open-loop system impulse residue magnitudes presented in Chapter III and 
those of Figure 4.9, one can see that this control law has decreased the contri- 
bution of the first elastic mode to the 9 t and 9 r outputs, but it has increased 
the contribution of the second elastic mode to these outputs. This is in opposi- 
tion to the desired result. This is an interesting result because this desired 
eigenvector selection procedure “worked” for the other modes. In this exam- 
ple, attempting to change a mode that was not contributing much to the rigid- 
body outputs, has yielded a control law that has caused this mode to contri- 
bute significantly to the 9 r output. 

As can be seen from Example 4.1, selecting the desired second elastic 
eigenvector to be the open-loop eigenvector and the desired second elastic 
eigenvalue to be close to its open-loop location yielded a much better result. 


Comments upon Response Amplitudes 

As noted earlier, when comparing the time responses obtained in Example 
4.1, the amplitudes of the responses decrease as the control weighting, p, 
decreases or different control laws result in different steady-state responses (or 
D.C. gains). The D.C. gain of a response due to a unit step can be determined 
in the following way. 

As was shown earlier, a system’s input-output transfer function relation- 
ship is given by 

i(s) = C[sl n - A]- 1 BUp(s) (3.15) 

For a stable system, the steady-state response due to a unit step input 
(u p (t) = one degree , u p (s) = 1/s) is given by 

= - [CA-'B] (4.21) 

For example, the steady-state angle of attack response due to a unit step 
pilot input is given by 

a ss ~ ~~ [ b ] — C a 

where 

c = row of C matrix associated with a 

b = column of B matrix associated with <5 e 

This equation shows that a commanded elevator deflection input of, for 
example, <5 e = (one degree) yields a steady-state angle of attack of C a g t degrees. 
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The larger the magnitude of C a the larger the effect S e has upon the steady- 
state angle of attack. Values obtained for a ss for the open-loop large flexible 
aircraft (Configuration Two) and several of the stable augmented systems 
presented ip examples in Chapters HI and IV are given in Table 4.10. 


Table 4.10 

Steady-State D.C. Gain 


System 


Open-loop 

-10.38 

Example 4.1,/?— 1 x 10 6 

-10.38 

Example 3.2 (DEA, six meas.) 

-6.893 

Example 4.1, p=l x 10~ 2 

-0.6293 

Example 4.1,p=l x 10~ 4 

-0.5388 

Example 4.1,p=l x 10~ 6 

-0.1460 


One may gain insight into how this occurs by calculating the steady-state 
value of elevator deflection due to unit step pilot input. Recall, the control law 
is given by 

11 = He + Sip = Gx. + ] 1 P (3.1e) 


where 

u c = feedback control law input 
u p = pilot' s input 

G = matrix of state feedback control law gains 
This equation can be rewritten in the following way 


K 


o 

X 

cT 

O 


M 

+ k 



^cv,xj v,8 t 


W 

lo. 


(4.22) 


where 
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— ( a > u f > ) £l > ?2 > £l > £2 > ^cv ) 

6 e = elevator deflection (scalar) 

8 ec = commanded elevator deflection (scalar) 

<5 ep = pilot' s commanded elevator deflection (scalar) 


The equation for <5 e is, for example, 


6. = G e>x X + G.*$. + 


e,S t °e 


(4.23) 


= «'. + G t,«A + 

Now, the elevator actuator dynamics are described by 

T T 


(4.24) 


Substituting Equation (4.23) into Equation (4.24) one obtains 


*. = (-7 + 


G 


r x. + 7(“'« + y 


(4.25) 


Taking the Laplace transform of this equation and finding the transfer function 
between 8 e and pilot’s input, <5 6p , ( with u' c = 0 ) yields 


_^e_ _ 1 

6 ep r s + (1 - G e<St ) 


(4.26) 


The steady-state value of elevator deflection due to a unit step in S e is 

1 


( 1 - g.a ) 


(4.27) 


This equation shows how the elevator deflection is reduced as the magni- 
tude of G 5<s increases. If it is desired to hold the D.C. gain of a steady-state 
response constant or meet some desired system output specification, one can 
simply multiply the pilot’s input by the constant ( 1 — G e ). 
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LQEA Conclusions 

LQEA is a technique for selecting quadratic weighting matrices that when 
used in an LQ synthesis will asymptotically place (n-m) finite eigenvalues and 
eigenvectors and m infinite eigenvalues and eigenvectors in the closed-loop sys- 
tem. With LQEA the designer has information concerning the asymptotic loca- 
tion of all of the eigenvalues and eigenvectors of the closed-loop system. As p 
tends toward zero in the LQ synthesis, (n-m) of the closed- loop system eigen- 
values tend toward desired locations. Their (n-m) associated eigenvectors tend 
toward the (n-m) system achievable eigenvectors. The remaining m closed-loop 
system eigenvalues tend toward infinity in first order Butterworth patterns. 
Their m associated eigenvectors tend toward the vectors E\°- 

Using LQEA, the designer determines which (n-m) of the system’s eigen- 
values are to be finite eigenvalues/eigenvectors. The remaining m are infinite 
eigenvalues/eigenvectors. One then chooses the desired eigenspace to obtain 
the desired synthesis objective. The desired eigenspace is used to determine 
the state weighting Q, and the control weighting matrix, R. These weighting 
matrices are used in an LQ synthesis with various values of control weighting 
to yield a family of control laws. 

The examples presented in this chapter illustrate the flexibility and poten- 
tial usefulness of this control synthesis technique. As with DEA, the proper 
choice of the desired eigenspace is an important part of the LQEA synthesis. 
The examples showed that the choice of the desired eigenspace can significantly 
effect the synthesis results. 

One potential disadvantage to this technique is that it yields a full state 
feedback control law. In many physical systems, the complete state can not be 
measured. Even in the cases where it can be, this type of control system may 
require too many sensors and may therefore be too complex and costly to 
implement. 

Also, in some cases, when it is necessary to use low values of control 
weighting to obtain the desired synthesis results, the LQEA technique may 
yield augmented systems with very high actuator bandwidth requirements. 
This could lead to problems with the practical application of this method. 

The next chapter will explore how some of the disadvantages of both the 
direct and LQ techniques can be overcome by using these techniques with state 
estimation. 
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CHAPTER V 

USING EIGENSPACE ASSIGNMENT TECHNIQUES 
WITH STATE ESTIMATION 


In chapter El, the Direct Eigenspace Assignment technique was presented. 
This technique was shown to be an attractive method for eigenspace assign- 
ment because it could be used to improve system dynamics and it yielded a 
relatively simple measurement feedback design. But as discussed, DEA had a 
number of disadvantages. In Chapter IV, it was shown that these disadvan- 
tages could be partially overcome by using the LQEA control synthesis tech- 
nique. But, overcoming these limitations comes at a cost of increased system 
complexity with full state feedback. 

This chapter explores how some of the disadvantages of both the direct 
and LQ techniques can be overcome by using state estimation. By using state 
estimation with the direct technique (full-state feedback) one has the freedom 
to place all of the modes of the system and one obtains the practicality of a 
measurement feedback control law. Using state estimation with the LQ tech- 
nique one can obtain a measurement feedback control law. 


Effect of State Estimation on System Dynamics 

How the state estimator dynamics affect the system transient responses, 
what method should be used to synthesize the state estimator, and how the 
state estimator affects the controller dynamics are three questions that will be 
addressed. The estimator dynamics could significantly affect the system’s tran- 
sient response and thereby nullify any effort expended improving plant dynam- 
ics. 

The effect of the state estimator dynamics on the system transient 
response can be examined in terms of their effect upon the augmented system 
transfer functions. This is done in the following way. 
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Given the system 


Ax + Bu + Ew x (system dynamics) 

(5.1a) 

X = Cx (system responses) 

(5.1b) 

= Mx + v z (system measurements) 

(5.1c) 

u,. = G£ (control law) 

(5.1d) 

H = ULc + lip (total control input) 

(5.1e) 


where 


x = estimated state 


Up = pilot' s input 

with x e R n , ue R l , £6 /? p , and xeR n . The vectors, w x and v z , are random 
processes reflecting external disturbances and sensor noises, respectively. 

Define the state estimator dynamics by 

x = Ax + Bu + F(z - Mx) (5.2) 

where F is a matrix of chosen estimator gains, and the state estimation error is 
defined by 

£ = X. - i (5.3) 

The system equations as functions of x and £ become 

£ = A& + Bn + Ew x 


= Ax + B(u c + u p ) + Ew x 
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= Ax. + BG£ + Blip + Ew x 

Therefore, 

x = (A + BG)k - BG£ + Bu p + Ew x (5.4) 

The error dynamics are obtained from Equation (5.1a) and (5.2), or 

£ = £~ £ 


= (Ax + Bu + Ew x ) - (Ax + Bu + F(z - Mi)) 


therefore, 

£ = (A - FM)£ - Fy z + Ew x 


(5.5) 


The augmented system then can be written 


X 


A + bg! -BG 


X 


B 


jc, 

— 

o |a-fm 


X 

• + 

0 

% + 


e!_o 

E | -F 



(5.6) 


One can now calculate the u p -to-x transfer functions after taking the Laplace 
transform of Equation (5.6). One obtains 


x(s) 


(sl n - (A+BG))! -BG 

-1 

B 

,i(s), 


o' ■ ; (si;-1aVmF) 


0 


i 

E* 0 


W 

E^ -F 




(5.7) 


(sl n - (A+BG))' 1 1 (sl n - (A+BG)) -1 BG(sI n - (A-FM)) -1 

0 ■ " (sl n -(A"FMr 

i 


B 

o' 
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E i 0 


W 

e;-f 


Zz. 


Since 


X = Ci(s) 


and assuming w x = 0.0, v z = 0.0; i(s) is given by 

i(s) = [sl n - (A+BG)] _1 BHp(s) 


(5.8) 


then, 

i(s) = C[sl n - (A + BG)]-'Ba p (s) (5.9) 


By referring to Equation (4.20) of chapter IV, one can see that the u. p -to-y 
transfer functions of this system with state estimation (Equation (5.9)) are the 
same as the u p -to-y transfer functions without state estimation. Therefore, as 
implemented here, the closed-loop system transient response due to a pilot’s 
input is independent of the state estimator dynamics. In this formulation, the 
estimator dynamics are uncontrollable from the pilot’s input, which is obvious 
from Equation (5.9). In the closed-loop system, the estimator poles (or the 
eigenvalues of (A - FM)), will always be canceled by zeros in each of the indi- 
vidual y 5 -to-Up transfer functions. A numerical example of this estimator 
pole/zero cancellation will be presented in Example 5.1 . Because of this 
independence, the state estimator can, in fact, be synthesized using any 
appropriate method. This yields design freedom to meet additional objectives 
such as robustness. For example, by using a Kalman filter as a state estimator 
and the method of Doyle and Stein [16], one could obtain enhanced closed-loop 
robustness properties. 

Unlike the closed- loop system transient response due to a pilot’s input, the 
choice of state estimator does affect the controller dynamics. The effect of 
state estimation on the controller dynamics can be explored by determining the 
controller transfer functions. These transfer functions .can be determined in the 
following way. Recall, the state estimator dynamics are given by 
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x = Ax + Bu + F(z - Mx) (5.2) 

= (A + BG - FM)x + Bu p + Fz 

= A con x + Bu p + Fz 

where A con = A + BG - FM. 

Taking the Laplace transform of this equation yields 

&(s) = (sl n “ A con r 1 Bii p (s) + (sl n - A con ) _1 Fi(s) (5.10) 

The control, u c , is related to i by the relation 

He = Gi 

Therefore, 

a c (s) = G(sl n - A co „r >Bu p ( S ) + G(sl„ - A co „r'Pi(s) (5.11) 
The total control is given by 

H(s) = u c (s) + u p (s) 

So, the controller transfer functions are given by 

ll(s) = [I„ + G(sl„ - A con )->Blu p (s) + G(sl„ - A con r‘Fi( 9 ) (5.12) 

From this equation one can see that the controller output is affected by 
both the measurements and the pilot’s inputs. The controller poles, the poles 
of A con , are directly dependent upon the state estimator. In the case of using a 
Kalman filter as a state estimator, the controller poles are dependent upon the 
choice of system noise intensity. A numerical example of this dependence of 
the controller dynamics upon the state estimator will be presented in Example 
5.1. 
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Control Synthesis Technique Summary 

In summary then, the solution method for this technique is a two part pro- 
cess: 

1) Use either the direct technique (full-state feedback) or the LQ technique to 
achieve the desired closed-loop eigenspace. 

2) Synthesize a state estimator using any appropriate method. 

An example will be presented in the next section to demonstrate the use of 
the LQ assignment technique with state estimation. 


LQEA with State Estimation Example 

How the LQEA synthesis technique can be used with state estimation will 
be demonstrated by an example. This example will extend the LQEA syn- 
thesis results obtained in Example 4.1 for p - 1.0 x 1Q~ 6 . A Kalman filter will 
be used as a state estimator. The model used is the same as that presented in 
chapter IV, Example 4.1 . The Kalman filter will be synthesized for the meas- 
urement vector 

1 T = ( ,&t f n 2 p » ) 

where 

0 t = total pitch angle (radians) 

0 t = total pitch rate (radians/sec) 

n Zp = plunge acceleration at the cockpit (g* s) 

n Zt = plunge acceleration at a tail station (g* s) 


These are the same four measurements considered in Chapter HI, Example 3.1. 


Example 5.1 

This example will start with the synthesis results of the LQEA example 
for p — 1.0 x 10~ 6 . 

The augmented system eigenvalues, the eigenvalues of (A + BG) were: 
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X sp = - 1.493 ± j 2.365 = ( w n = 2.797 (rad/s) , 0.534 ) 

X ph = - 0.0075 ± j 0.0679 = ( w n = 0.068 (rad/s) , 0.110 ) 

X ?i = - 1.807 ± j 8.608 = ( w n — 8.796 (rad/s) , ?= 0.205 ) 

X & = - 0.5351 ± j 21.35 = ( u n = 21.36 (rad/s) , 0.025 ) 

Xgl — 500.5 

X cv = - 999.9 

The eigenvectors of (A + BG) can be found in Chapter IV, Table 4.7. 

The Kalman estimator will be synthesized for a fictitious state noise inten- 
sity (actually a design variable selected to yield an acceptable design) of 

W = I 

and a measurement noise intensity (also, in effect, a design parameter) of 

V z = «I m 

where k is a positive scalar quantity. 

The estimator will be synthesized for two values of k (1.0 x 10 6 and 1.0 x 10 -6 ) 
to demonstrate how the choice of state estimator can significantly effect the 
controller dynamics. 


Control synthesis for k = 1.0 x 10 6 . The filter gain matrix and the augmented 
system matrices for this value of p and k can be found in Appendix A.5. The 
system estimator eigenvalues, the eigenvalues of (A — FM) are: 

X 12 = “ 0.4594 ± j 21.35 = ( w n = 21.36 (rad/s) , f= 0.022 ) 
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X 34 = ~ 0.7268 db j 8.758 = ( w n = 8.789 (rad/s) , 0.083 ) 

X 5 6 = - 1.349 ± j 2.193 = ( w n = 2.574 (rad/s) , ?= 0.524 ) 
X 78 = - 0.0042 ± j 0.0525 = ( w n = 0.053 (rad/s) , 0.079 ) 

Xg “ — 10.0 Xio - ■ 9.0 


The eigenvectors of the complete augmented system are given in Appendix A.5. 
Note that the eigenvectors associated with the poles of (A + BG) are the same 
as before. 

The poles of the controller, the poles of A con are: 

X 12 = -0.538 ± j 21.35 

X 34 = -1.809 ± j 8.610 

X 5 6 = -1.495 ± j 2.366 

X 7 8 = -0.0052 ± j 0.0679 

X 9 = -999.9 X 10 = -500.5 

To demonstrate the estimator pole/zero cancellation in the y-to-u p transfer 
functions, the 0 r -to-6 e transfer function is determined (recall, in this model 
u p T = ( <5 e , 0.0 )). The poles of this transfer function are the eigenvalues of 
(A + BG) and (A — FM). The zero’s of the 0 r -to-<5 e transfer function were cal- 
culated using the method of Sandburg and So [14]. This transfer function is 
given in Figure 5.1. As can be seen, the poles of (A - FM) are a subset of the 
transfer function zeros and will therefore cancel in the transfer function. 
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Figure 5.1 
d T /6\ Transfer F u 
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The output time responses due to a unit step in u p , with zero initial condi- 
tions, are shown in Figure 5.2. As is expected, these time responses are identi- 
cal to those obtained in Example 4.1 for p — 1.0 x 10 6 (See Figure 4.8). 

Control synthesis for k = 1.0 x 10 -6 . The filter gain matrix and the augmented 
system matrices for this value of p and k can be found in Appendix A.5. The 
system estimator eigenvalues, the eigenvalues of (A — FM) are: 

X 12 = - 2.244 ± j 2.035 = ( u n = 3.030 (rad/s) , f= 0.741 ) 

X 3 = - 1.106 x 10 5 X 4 = - 7.326 x 10 4 

X 5 = - 1705.7 X 6 = - 1606.1 

X 7 = - 96.95 X 8 = - 9.974 

X fl = — 6.209 X 10 = - 1.566 

The eigenvectors of the complete augmented system are given in Appendix A.5. 

The poles of the controller are: 

X 12 = -0.0017 ± j 75.98 
X 3 = -1.106 x 10 5 X 4 = -7.35 x 10 4 
X 5 = -781.5 X 6 = -411.2 
X 7 = -54.35 X 8 = -8.249 
X 9 = -2.497 X 10 = 0.0728 



GAMMA (RAD) CX10 ~ 2 ) „ U (FT/S) 
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Figure 5.2 

Output Time Responses - Example 5.1 


THETfl-R DOT (RflD/S) (XIO ~ 2 ) THETA-R (RAD) (XIO " 2 ) 
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Figure 5.2, continued 




THETR-T DOT ( RRD/S ) (XIO " 2 ) THETfl-T ( RRD ) (XIO ~ 2 ) 
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Figure 5.2, continued 
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TIME (SI 


Figure 5.2, concluded 
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ORIGINAL PAGE 13 
OE POOR QUALITY 


The output time responses due to a unit step in y. p , with zero initial condi- 
tions, are identical to those of Figure 5.2. 


Comments on Example 5.1 

This example illustrates how the LQEA control synthesis technique can be 
used with state estimation. It is significant that the regulator gain synthesis 
and state estimator synthesis are independent. Therefore, the impulse residues 
and time responses for each case here are the same as the impulse residues and 
time responses for the p — 1.0 x 1CT 6 LQEA synthesis. 

The Kalman filter synthesis for k — 1.0 x 10 6 is a synthesis for a large 
measurement noise intensity. This synthesis produces small filter gains and an 
estimator with a small bandwidth (poles of (A — FM)). For this synthesis, the 
controller poles are stable. Because of the small filter gains, the poles of the 
controller are very close to those of (A + BG). 

The Kalman filter synthesis for k — 1.0 x 10 -6 is a synthesis for a srn 7 
measurement noise intensity. Therefore, this synthesis produces large filter 
gains and an estimator with a large bandwidth. As can be seen from the con- 
troller poles, the controller is unstable for this synthesis. From a practical 
point of view this is very undesirable. 

Conclusions 

This chapter explored the advantages of using the direct and LQ eigen- 
space assignment methods with state estimation. By using either the direct or 
LQ methods with state estimation one can control the placement of all of the 
poles of the closed-loop system and obtain a perhaps more practical measure- 
ment feedback control law. 

As was shown in the development and in the example, the closed- loop sys- 
tem transient response due to a pilot’s input is independent of the state estima- 
tor dynamics. This is because in this formulation, the estimator dynamics are 
uncontrollable from the pilot’s input. Because of this independence, the state 
estimator can be synthesized using any appropriate method. 

Although the choice of state estimator does not effect the a. p -to-y: dynam- 
ics, it does effect the controller dynamics. As was shown in the examples, some 
choices of state estimator can even lead to an unstable controller. Therefore, 
the controller dynamics must be considered when choosing and synthesizing the 
state estimator. 
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CHAPTER VI 
CONCLUSIONS 


This report has presented two eigenspace assignment techniques: Direct 
Eigenspace Assignment and Linear Quadratic Eigenspace Assignment; and has 
explored the advantages of using these techniques with state estimation. These 
techniques have been applied in several cases to the problem of control law syn- 
thesis for a flexible aircraft. 

Direct Eigenspace Assignment (DEA), presented in Chapter HI, is a control 
synthesis technique for directly determining measurement feedback control 
gains that will yield an achievable eigenspace in the closed-loop. For an 
observable, controllable system that has n states, m controls, and l measure- 
ments one can determine a gain matrix that will place / eigenvalues to desired 
locations and their / associated eigenvectors as close as possible in a least 
squares sense to desired eigenvectors (this assumes / > m). This technique was 
shown to have two major disadvantages. First, there is no systematic way to 
trade eigenspace assignment goals versus gain magnitudes. Second, the (n-f) 
unspecified system poles can move to unpredictable locations, sometimes even 
destabilizing the system. 

These disadvantages were partially overcome in the Linear Quadratic 
Eigenspace Assignment (LQEA) control synthesis technique. LQEA is a syn- 
thesis technique for selecting quadratic weighting matrices that when used in 
an LQ synthesis will asymptotically place (n-m) finite eigenvalues/eigenvectors 
and m infinite eigenvalues/eigenvectors in the closed-loop system. It yields a 
guaranteed stable closed-loop system. The designer has information about the 
asymptotic behavior of all of the modes of the system and has the ability to 
trade eigenspace assignment goals versus gain magnitudes. But, using LQEA 
comes at the cost of increased system complexity due to the necessity of full 
state feedback. This is a problem because many times the complete state can 
not be measured. Also, LQEA may yield augmented systems with high actua- 
tor bandwidth requirements. 

Chapter V explored the advantages of using the direct and LQ eigenspace 
assignment techniques with state estimation. By using state estimation with 
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either technique the designer has control over the placement of all of the poles 
of the closed-loop system and can obtain the practicality of a measurement 
feedback control law. 

As was shown in Chapter V, by using the proper formulation the closed- 

loop system transient response due to a pilot’s input is independent of the esti- 

mator dynamics. Therefore, the state estimator can be synthesized using any 
appropriate method. But, the choice of state estimator does affect the con- 
troller dynamics. Unlike the closed-loop system, the controller is not 
guaranteed to be stable. The controller poles are directly dependent upon the 
state estimator. Some choices of state estimator can result in an unstable con- 
troller. This is very undesirable from a practical point of view. Therefore, the 
controller dynamics must be considered when choosing the estimator. 

There are many areas for future work. Some of these include: 

1) The art of choosing the desired eigenspace. As was pointed out in the 

examples, the choice of the desired eigenspace is very important and can 
significantly effect the synthesis results. Much work can be done in the 
area of exploring what desired eigenspace is needed to achieve a certain 
desired dynamics. For example, what desired eigenspaces would yield 
good handling qualities? 

2) Extending DEA to be able to trade gain magnitudes versus eigenspace 
assignment. DEA would be much more useful if one knew how a change 
in an element of the desired eigenspace would affect an individual gain 
magnitude or the location of one of the unplaced poles. 

3) Exploring how the phase relationships between desired eigenvector ele- 
ments would effect the synthesis results. Choosing proper phase relation- 
ships could make the achievable eigenspace closer to the desired eigen- 
space. 

4) Determining under what conditions an LQG control synthesis yields an 
unstable controller. 

5) Investigating how the choice of measurements (sensor locations) effect the 
synthesis results. 

6) Determining how well each of these synthesis methods rejects distur- 
bances. 

7) Examining the sensitivity of each synthesis technique due to variations in 
the model parameters. 
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Appendix A.l 
DEA Development 


Case of Control Feed-Throuffh in the Measurements 
Given the system 

x = Ax. + Bu (system dynamics) 

(A.l.la) 

y = Cx + Du (system responses) 

(A.l. lb) 

1. — Mi + Nn (system measurements) 

(A. 1.1c) 

u c = Gz (feedback control law) 

(A.l. Id) 

u = u c + Up (total control law) 

(A. l.le) 

Up = pilot' s input 

where xeR n , u eR m , y eR p , zeR 1 . 

Solving for u c as a function of x yields 

He = G z 

(A.l. If) 

= G[Mx + Nu c ] 


= GMx + GNu c 


u c - GNu c = GMx 


He = |I m - GN]-‘GMx 

(A. 1.2) 

The closed-loop system is given by 

i = Ay + Bn 


= Ax + B[I m - GNp’GMx + Bu p 
x = (A + B[I m - GN] _1 GM)x + Bu p 

(A. 1.3) 
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Thr. spectral decomposition of the closed-loop system is given by 

(A + B[I m - GNr'GMjj/j = X^ i = 1, n (A. 1.4) 

Let W; be defined by 

Wj = [I m - GNpGMt* (A.1.5) 

Substituting Equation (A.1.5) into Equation (A. 1.4) and solving for E\ one 
obtains 

Ei = (Ail n - A)“‘Bwi (A. 1.6) 

Equation (A. 1.6) is an expression for the achievable eigenvectors of the system 
given Xj and Wj. Note that this is the same result as obtained in the N = [0] 
case. An expression for w ; can be obtained by substituting (A. 1.6) into the cost 
function 

Ji - \ - I'd) (A. 1.7) 

where 

E a ( = i th achievable eigenvector 

Ed { — i th desired eigenvector 

Q d . = i th weighting on eigenvector elements. 

Taking the gradient of J; with respect to w i( setting this equal to zero, and 
solving for Wj yields 

*i T = "d, T Q - Ar 1 B[B T (X d J„ - A)- T Q d ,(X di I n - A)B]-> (A.1.8) 

Concatenate the individual w;’s and i/j’s columnwise to form W and V 
respectively. Solve for the gain matrix, G, using Equation (A.1.5) 

W = [I ro - GN] _1 GMV (A.1.5) 

[I m - GN]W = GMV 
W = GMV + GNW 


G = W[MV + NW ]” 1 


(A. 1.9) 
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Complex Eigenvalues and Eigenvectors 
System 

x = Ax + Bu (system dynamics) 
z = Mx + Nu (systemmeasurements) 
H = + Up (total control input) 

U c = Gi (feedback control law) 


F or this case the cost function is given by 

J i = J K ~ ~ I'd) (A. 1.10) 

where superscript * denotes complex conjugate transpose. Again, the spectral 
decomposition of the closed-loop system is given by Equation (A.1.4) 

(A + B[I m - GNp'GMlj-i = X,n, (A.1.4) 

Let 

El = £r. + m, (A.1.11) 

Xi = <Ti + jo/j (A.1.12) 

Equation (A.1.4) can then be written as 

(A + B[I m - GN^GMXifc, + in,,) = (<T, + jUiXnk, + inj (A.1.13) 

and its complex conjugate 

(A + B[I m - GNl-'GMXih, - in,,) = (<t, - juOlnh, - jn,,) (A.1.14) 

Multipling out either of these equations one obtains 

(A + B[I m - GN]-'GM)n E , - i(A + B[I m - GN]-‘GM)n,, 

= o',!*, - inR", - in,,*, - "in,, (Ai.15) 


Separating the real and imaginary parts one obtains 
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(A + B[I m - GNpGM)^. = ^ Rj - ^ (A.1.16R) 

(A + B[I m - GNp’GM)^. = i/R.Wi + mp\ (A. 1.161) 

In these equations let 

SEr, = (I™ - GNr'GM £Rl (A.1.17R) 

ffil, = (t m - GNr'GM £li (A.1.17I) 


Substituting Equations (A.1.17R) and (A. 1.171) into Equations (A.1.16R) and 
(A. 1.161) yields 

Al/R. + Bw r . = (Tji/R. “ 0J&. 

A Ei, + Bwj. = i/R.Wi + <?■&. 

In matrix form this becomes 


b|o 


w Ri 


((Tjl n A) , 


o ;B 


as. 


I (^n " A ). 



Bw- X — M(\ j)^j 


(A. 1.18) 


The expression for the achievable eigenvectors of the system is given by 


Mi = AT^X;)!^ 


(A. 1.19) 


where 


— Ri 

w zz 

Wr, 


) _ 1 

^Ii 


The equation for w x is determined in the same way as before. For this case the 
equation becomes 
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»i T = B T Ar T (x)Q dl Ar'(X)B ]-' 

where 





< 0 

= ' 


> Qd;~ 

~°~\S 


The gain matrix is solved for using Equations (A.1.17R) and (A.1.17I). 

W = (I m - GN) -1 GMY 


where 


G = W[MV + NW ]" 1 

w = [w Rl | w It | w Rz | w Ia | ... ] 
V = [v Rl | V lt | Vr 2 | V l2 | ... ] 


(A. 1.20) 
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Appendix A.2 
Open-loop System. Data 
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Figure A.2.1, continued 
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Appendix A.4 
LQEA Examples Data 
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Figure A.4. l 

Open-loop System Matrices - Modified Configuration Two 
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Figure A.4.2 

Spectral Decomposition - Modified Configuration Two 
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IMPULSE RESIDUE MATRIX ( 7 BY 10) 

OUTPUT FORM: ( MAGNITUDE , PHASE (DEG 
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Impulse Residues - Modified Configuration Two 
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Figure A.4.5 
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Figure A.4.7 

Closed-loop System Matrices - Example 4.1, 
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Figure A.4.8 

Spectral Decomposition - Example 4. 
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Figure A.4.8, concluded 
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Figure A.4.9 

Impulse Residues - Example 4.1, p — 10 
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Figure A.4.10 

Closed-loop System Matrices - Example 4,1, 
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Figure A.4.12 

Impulse Residues - Example 4.1, p — 10' 
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Figure A.4.13 

Closed-loop System Matrices - Example 4.1 
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Figure A. 4. 14 

Spectral Decomposition - Example 4.1 
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Figure A.4.14, concluded 
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Figure A.4.15 

Impulse Residues - Example 4.1, p = 10 
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Appendix A.5 
Chapter V Examples Data 
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Figure A.5.1 

Augmented System Matrices - Example 5.1, 
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Figure A.5.1, continued 
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Figure A.5.1, concluded 
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Figure A.5.2 
Spectral Decomposition - Exa 



THE EE MATRIX ( 20 BY 20 


209 


ORIGINAL PAGE IS 
OF POOR QUALITY 


If) rx r- » /-v ^ /-» /— x /-v ^ rv r\ 

nmnnnjinHnno 
oooooooooo 
1 I I 1 I I I I I 4 

LJUULdUliJLdLiJbJLjJ 

^rcn^iooooTtDxrruo 

coinooioj-^-rucoin^ 

nCOHHtHOOntlOHOOOOOOOOOO 
II I I I I I 

co co in ^ ^ t? h ru *h 
oooooooooo 
l I I l l l l I I t 

LUliJLLHjJUiLJLUUiLiJIjJ 

Trcnto^-totoaicooLn 

lo^ranoojHinin^ 

ntnHNHt(DCOn(\Joooooooooo 
I I I I I I I 


W W W W W ' 


' W Vj» V/ V 


- w w w w v-/ w 


rv ^ / 


«-N /— v rx y — v rs rx y — v rx y 


cociJtosrtoaicuoco* 
oooooooooo 
l I I i l t I 4 I I 
LJLdUJUJLdiiJUJhJbJljJ 

cooocdcooididodcd-h 

-HCo-*\rcnintD^ina) 


HH(lJKniH{\|HU)HOOOOOOOOOO 

till 111 


^rrucxi^-coajnj— «roru 

OOOOOOOOOO 

I 1 1 I I I I I I I 

UJLJULjJliJLdlULJLiJliJ 
in^rruN-N--H(oo-njro 
incoLncn^in-HoiDOo 


inHtutOHinmnaj^oooooooooo 

i i i 

w w w vy w w w w w w w w w w w' w w w w 


fQ /A rx rx /x rx rx /-v r\ rx rx y 

naju 34 nainjon-< 

OOOOOOOOOO 
I I I I I I I 4 I I 
IxlLULJLJbJUJLJLiJLiJljJ 

roootofooamcocn-H 

^HCO-«^cnLntD-untD 


— I— tCUN-OJ— iCU-hCO-hoooooooooo 

I I I 


TrnjcD^majru-Hforu 
oooooooooo 
i i i i i i i i i i 
UJIxJUJUJLiiliJUJliJUJtii 
inxrajr^f^-Hro^ojn 
iDCDincn-an-Hoinco 


ir)Hcuc0Hinu)najtoooooooooo 


i i 

X-/ W W WWW XW W W W W W W W W x 


ru ' 


x rx rx rx rx rx rx / 


oooooooooooooooooooo 


TrfucoLnxrLnnjnj— <nj 
oooooooooo 
i i i I i i i i i i 
UJUJUJUJLUUJUJLiJLiJliJ 
COCOCDCDN-CnN-LnofO 
-HCOCOh-COCDCOCOCO-* 


^OJH'tHtDOinMDoooooooooo 

III I 111 

W W >W W- W WWWWWWWW «W WWW W' w x_/ 

I r\ /-x /~x rx /-x /x y-v rx rx rx rx r\ rx rx rx yx rx /x yx 


oooooooooooooooooooo 


^cncncoinN-OJ^-i-H 
oooooooooo 
I I I I 1 I I I I 4 
IxILUUJLlILlILlIUJUJUJIjJ 
o^Hin-HCocncocntntD 
oionjoN-nh-ncDin 

'-tmrjoo'tfTQTfn— i-hoooooooooo 

i I i 


) yx y-x rx rx y-X rx rx rx rx rx / 

iruru-»-Hronin\rn-« 
oooooooooo 
l l l i i i ( i t i 

LULULJLiJliJLiJliJLULjJLJ 
4 ^imnncDtnootj* 
oironcn— iooocq^it) 


(nnainwtajmnoooooooooo 

ii r i i i t 


ojoj— 

oooooooooo 
I I 1 I 1 I ! I I 1 
IxJUJUJLULJUJLlJLiJLiJIjlI 
rrco^-cn^-^cnooo-H 

co ro (v^^tnN-CDcoco 


NOJHtHnain^fuoooooooooo 

i I 


m - 


> yx yx rx / 


ojoj-H—froDin^m-n 
oooooooooo 
l l I l I i I I i l 
UJUJljJLdLiJLiJUJUJUJLiJ 
rr^uncnrocDcnooT}- 

ru ro co tn ^ co o co o- in 


(nnajnt^tnjinnoooooooooo 

i I i 


runj— »-i^roo-o-nj-H 
oooooooooo 
1 I I I I I 1 1 1 1 

IxJJjJUJUJLlJIJLlJIjJLiJLlJ 

o-cooroo-o-cnoco-H 

coror^^TFcnrvtococo 


M\|H^H{ 00 JnH(\|OOOOOOOOOO 


I 


QQ yx rx rx yx rx rx yx rx rx rx/ 

^^Hrrj^cnojfo^njo 
oooooooooo 
t i I i i I i f I 4 
UiLULiJUJUJUJLJlUUJUJ 
— iN-ino-omimncnru 
-HCDCDo\rruooojaj 


/ w w -w W W W> W w w 


, rx rx rx 


' w w w 


I y-x rx rx f\. j 


OJfOlDHOiHfOHtHoooooooooo 

111(1 

-h j; ^ j; xj co co ru ru 
oooooooooo 
I I I I I I I I 1 I 
(jJIxJIxILiJLiJUJLlIIjJLJLJ 
*H 4 inr)C 0 Tj-njinon 
cncnaiforucno-Hncn 

HinnaiajKUJ^Hnoooooooooo 
t it i 

w x_y www w w w w w w w •w w w x-y ■w w w w 


oooooooooo 
I I 1 I 1 I I I i 4 
LJLdUJLiJLiJLiJLLlbJLdLiJ 
^Mn^oimnmmaj 
HoitDo^ojooaJOJ 


. y-x rx y-x rx rx rx rx rx rx rx rx 


(UfoinHajHCOHtHoooooooooo 
i i i i i 


oooooooooo 
l l I i l i l 1 I l 
UJUJIxIbJLdljJLiJLiJLiJLiJ 
-t^rinrncoTfnjinorQ 
cnmojnajtno— (men 


HinD(\j(UMO^Hnoooooooooo 
i ii i 


f n rx rx rx rx rx rx rx rx rx rx rx 

cocoLTMoaiin-tcocoo 
oooooooooo 
I I I I I I I I 1 4 
LUIxJLjJLlJIiJIjJUJLjJLlJLiJ 

itococnto'o-njo 

coinoaico^rrucDLD-H 


rx rx rx rx rx rx , 


ntOHHtHOOfOtOHOOOOOOOOOO 


I 


nmn 4 ^ 4 w^-njH 
OOOOOOOOOO 
I f I I I I I I I I 
UJUJUJbJLiJUJljJtiJUJliJ 
^CDCDxrCDLDOJCOOLJl 

ui^rcocoonj-<tnino- 


nmHKHttflCQnojoooooooooo 

i i i i i i i 


^njn^rmcDr^oocno. 


<airo^riDLDN.co(Do 

4 ■— I r — 4 ■— I - — 4 <H *— 4 >H OJ 


jnjro^rinuDivcocno—tnjn^rinuDr^cocno 

hhhhhhhhhh[\J 


Figure A.5.2, continued 
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Figure A..5.2, concluded 



S0-3000000T * = M>X ++++++++++ 

S0-300000QT* = QHd 803 dOO - !********** 


211 


I ^ 

I 

i 


> 

PQ 


I o 


I w 

I 

I X 

I *H 

I £K 
I h- 
I <X 
I 21 


I U- 


UJ 

X 


TT cu 

cu 

cu 


•H 

ru 

o 

o 

o 

o 

o 

o 

+ 

+ 

+ 

+ 

+ 

+ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

CO 

•0- 

cu 

o 

LD 

cu 


07 

o 

OJ 

N. 

in 

Tf“ 

cn 

CD 

07 

O 

in 

o 

TT 

CD 


CO 

07 


m 

*■ 

• 

• 

9 

iv 

CD 

•nT 

N- 

1 

ru 

rt 

CO cu 

H 

cu 

•H 

CU 

CU 

o 

O 

o 

o 

o 

O 

+ 

+ 

+ 

+ 

+ 

+ 

UJ 

Id 

UJ 

Ld 

Ld 

Ld 

07 

CD 

cu 

cu 

CU 

CO 

o 


07 


CO 

H 

07 

CD 

in 

in 

IV 

LO 

N- 

rv 

»-* 

CO 

m 

07 

n! 

H 

i 

i 



i 

CU -H 

cu 

cu 

cu 

cu 

cu 

o 

o 

o 

o 

o 

o 

+ 

+ 

+ 

+ 

+ 

+ 

Ld 

u 

UJ 

Id 

UJ 

Id 

CO 

CO 

CD 

07 

fv 

LD 

o- 

in 

\r 

O' 


07 

o- 

cu 

07 

CO 

07 

CD 

CO 

CO 

CO 

07 

CD 

CO 



« 




OD 

cu 

to 

rv 

in 

CU 

sH H 

cu 

cu 

cu 

cu 

CU 

o 

O 

o 

o 

o 

O 

+ 

+ 

+ 

+ 

+ 

+ 

id 

id 

UJ 

Ld 

Ld 

Id 

o 

CO 

CO 

CO 


CO 

07 

07 

CD 

iv 

CO 

H 

CO 

O 

in 

O 

r- 

00 

H 

in 

cu 

CU 

O 

in 

• 

• 

0 

• 

b 

• 

cn 

cu 

CO 

LD 

N- 

1 

cu 

l 


cu 

CO 

o- 

in 

L0 


cucoru-t 

0 o o o 
+ + + I 
Ld Ld Id Id 

01 co cn 
CO o N- CD 
Kinmo 
in in to 


f oj m* 


+ + + + 


o co h- in 


i ai — » ai 
t 

main 

> o o o 
+ + + 

1 UJ UJ UJ 

> OJ LD 07 

j co otn 

' O CO h 
' N- CO CU 

I « • » 

* OJ 
I I 


i o 



i <r 

OJ to — 1 o 

1 21 

o o o o 

1 

+ + + + 

1 

Id Id UJ Id 

1 

oo N- co ru 

I 

o in cn co 

1 

CO CD CO CU 

l 

N. ’•hi O CO 

1 

• • • • 

i a 

CO-i-tf* TJ- 

1 3 

1 1 

i <r 


1 <E 

N. CO 07 o 

1 Ld 

•-I 

l X 


l »“ 


07 •H •H 

hhCUh 

oo 

o o o o 

1 + 

+ + + + 

Id UJ 

Ld Ld id Ld 

o o 

CU LD IV 07 

CO 

N O 07 

070 

00 CU o OJ 

coin 
• • * • 

07 •h O CD 


WHOoooinaJinaioooooooooo 
i I III 


00 CO *h 

0 o 

1 I 
UJ UJ 

in oo 
in^r 

H H 

rt o 


OHHOtU 
o o o o o 
+ l I + + 
UJ Id UJ LU UJ 
OlDOKCn 

otoincuco 
ocncoinrv- 
o cn cn o 


^tuooo— IU3C00I-HOOOOOOOOOO 

It It 

N- CO -h 

o hhooi 

O o 

o o o o o 

t 1 

+ 1 1 + + 

Ld Id 

Ld Id Id Id Ld 

ru o- 

O O CD CU IV 

07 tv 

O CU U7 CD CU 

LD 

O 0107 07 U7 

O 00 

O LD LH CU -h 

coaioorioijjHcocnDooooooooo 

i 1 1 1 

CD *h O 

h cu o ru 

o o 

o o o o 

t + 

+ + + + 

Id Id 

Id Id td Ld 

cu 

o ^ tj- 

ld in 

CD in CD 07 

CO 07 

IV LD CO o 

(0 07 

•h LO «h 03 

•Hinoooo^-txj-anvoooooooooo 
It t I 

inno 

HOCUC7 

o o 

o o o o 

1 + 

+ + + + 

Ld Ld 

UJ UJ Ld Ld 

o- ro 

o in IV (O 

CD [v 

o o in in 

CD N- 

07 in CON- 

CO 0) 

CD — 1 o CU 


<o 

l 

o 


C0CJQooooivD7*Hfvoooooooooo 

1 f 

o- cu 

O *H 

o 

o o 

1 

+ 1 

UJ 

Ld Ld 

o 

CO 00 

cu 

ru co 

IV 

co 

CO 

coco 

o o CO o o o 
1 

OOCU07OOOOOOOOOO 

I 

co -h o ru 

ajHoaj 

o o o 

O O O o 


I + 1 
UJU1UI 
K OJ O 
CD CO O 
^ O 

o — « in 


i t + + 
Ld UJ UJ UJ 
nj-tcnro 
mcnoco 
07 rv iv m- 
lO OJ O 


cu t 

CD CU 

ru o o o 

1 III 

> 1 

PQ t 

CUHO 

o 

t 

O o 

o 

O ! 

1 + 

+ 

ru i 

Ld Id 

Ld 

1 

o o 

o 

W 1 

OCQ 

o 

1 

Tf 

o 

X 1 

NT o 

o 

1— 1 1 

• • 

« • • • 

Ck: 1 

mcuoHoo 


i i 

OO-tCU 
o o o o 
+ + + + 
Ld Id UJ Ld 
*h LD CO CO 
07 ro co 
*-f CO CD 
tn 07 N- cu 


l 

HOOOJ 

o o o 
+ + l 
UJ UJ UJ 
oo^r 
00 cn rt 
o in n 
CXJO LD 


i I 

-* CU to 
o o o o 
+ + + + 
Ld Ul UJ UJ 
-h LD h- CO 
07 CO O- 07 

xr -« cur^- 

07 LD O CU 


— iJvOJOOO—JfU— f’tf'OOOOOOOOOO 

(If III 


tnjcoo'intorvcocno-iaicosrinioivcocno 
H H H •— I H H •— I •— f *H «— * CU 


Figure A. 5. 3 

Augmented System Matrices - Example 5.1, 



-4.4000E-03 


212 


03 


o 

OJ 

cu 

cu 

ru 

ru 

cu 

CU 

ru 

ro 

CU 

—1 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

Ld 

Id 

Id 

id 

Ld 

Id 

Ld 

Ld 

Id 

Ld 

Id 

Ld 

N- 

cn 

in 

Tf 

in 

CD 

o- 

<v 

N- 

O 

TT 

in 

ru 

co 

cn 

UD 

rv 

h- 

ro 

in 

(O 

•H 

H 

M- 

in 

fv 

rv 

N- 

in 

in 

tn 

IV 

03 

03 

o 

w 

03 

• a 

o 

• 

xr 

• 

O 

• 

o 

« 

to 

* 

» 

o 

• 

00 

IV 

o 

« 

CO 

• 


03 


in 


CO 


Of 


ooooooooaiHW-iMDUjHconj^H 
i ill i it 

orururucurururucutocu-H 
i oooooooooooo 

++++++++++++ 
LdLdLdldLULdldLdldUJLdLd 
ruiv— itinviv^Lnruvcoo 
uoru-n^-ocncninocicncoco 
cmnnincoincDo-to^co 
OJ^o—iCD\rcocD^o^in 

OOOOOOOOtOtnN-NHNtD^HHHH 

i III 

oajTf^^n^^inmajn 
oooooooooooo 
++++++++++++ 
LdLdLdldLdLdldLdldLdLdLd 
^r^cnujN-^-njroinnncn 
cocnaDin^— JOiivcooDcnco 
rooLnuD^LDtoiocncocoo- 
— iCOCQN.COOiKf^tDON .00 

oooooooocuivnjtocoaDcocnajcotD-H 

till ti 

ruro^^Tf-HTj-^ininnaj 
oooooooooooo 
+• + + + + + 4 - + + + + + 
LdldLdldLdldldldldldLdld 
rvcocotuiDcoinroivtnTr*--* 
inincnrocococu-nN-cootn 
coi^tDODcocn^rinooo^r 
onjnjoocxjnjcn-ioonco 

oooooopoH|CnVair)HnHHU)co 
ti i i i t i 

o^-icuruajairunjco— «o 
* oooooooooooo 

+ I + + + + ++ + + + + 
LdLdLdLdLdLdLdUJLdLdLdld 
cooooooruroo-cocoivcoru 
oirocncnKh-n-Joincnn 
co + moinoMomoooai 
coco^incutuoinrv— ioro 

oooooooonjcomnjncDNajcoH^^ 

i lilt ii 

ocuo-'tf-rucoroioo-roronj 

OOOOOOOOOOOO 
++++++++++++ 
IdldidldLdLdldldldldldUJ 
(Dtnooo + HHtnnMn 
ocotocorucuTj-ivoorvtu-i 
k^ninmofontoMno 

OH + fUHHNO^OJOJO 

oooooooonajaiHinHtflHHncon 
ii till 

-Htucocococutucocoronjo 
oooooooooooo 
++++++++++++ 
LdLdLdldLdLdLdLdLdldldld 
COCOCOCO-HtDOtOLDCOCOO 
tDHMU^NOKOJntDtO 
-HintnmajrHfoin^inajtD 
N.nj^-H^incoatCDN-^r^r 

ooooooooHMvncUTrnjntDMnN- 
i i ii i 1 it 


o o o o o o 
++++++ 
ill Id Ld UJ UJ UJ 
N- CO (O N- LD LD 

cn cn cn — « to 
rur^ co inru in 
ociunnoo 


u- to 

o o o o 
+ + + + 
Id Ld Ld Ld 

in -h cn cn 

— * *h U3 CU 

in in rrn 

CO CO 00 o 


o- ru 
o o 
+ + 
Ld Ld 
LQ —t 

O tn 
jv co 
co 


ooooooooHtinto^nHnHHHin 

ill i i i i i i 


0 

1 

Ld 

o 

H 

CO 

in 


OhOCO 
o o o o 
+ I + + 
Ld Ld Ld Ld 
CU OJ CD in 
to in in ru 
to ru ru o 
+ ono 


(Oooooruann-Hoooooooooo 
i t 


o 

ru 


o oj to 
o o o 
+ + + 
Ld Ld Ld 
CO K CO 

in^rtn 

aicu-H 

ro co co 


rucoo 
o o o 
+ + + 
Ld Ld Ld 

to cn ru 
cn ru iv 

o- to 
- — • rxi — * 


euro 
o o 
+ + 
Ld Ld 
CD O 
tT 

o cn 

0J to 


to co ru 
o o o o 
+ + + + 
Ld Ld LJ Ld 

cn ru cn o- 

LD 00 
L0 CD IV O 

in co iv ro 


cn 


ooo oooooinoiHfijHaicnH^Tfnin 
ii ii i it 

ru-H^-o-o-coroo-^ro-co-H 
oooooooooooo 
++++++++++++ 
IdLULdLdldLdLdLdLdLdLdLd 
N-CnC 000 tDO|vnj{VfVN.C£l 
xt'cncoN-ruincncomtotoo 
<~icuoN-OQcnrurvrv(OU3cx3 
cntocn^-LDLoncDcnrurvcD 

ooo ooooo + njojnj-Hajai^alincoK 
til ii ill 


ncl 


^cucoTrinu)ivcocno»-iruro^rinu3ivcocno 

H H H H H *H *— I •— I H i — I QJ 


ru ro TMna)ivcotno^(uco-d-inu3N.cocno 

H H H H H H H H H H f\J 


Figure A.5.3, continue 
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Figure A.5.4, continued 
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Figure A.5.4, concluded 
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determining constant measurement feedback gains that will yield a closed-loop 
system eigenspace "close" to a desired eigenspace. The second technique is a 
method for selecting quadratic weighting matrices in a linear quadratic control 
synthesis that will asymptotically yield the closed-loop achievable eigenspace. 

Finally, the possibility of using either of these techniques with state estimation 
is explored. 

Application of the methods to synthesize integrated flight-control and struc- 
tural-moae-control laws for a large flexible aircraft is demonstrated and results 
discussed. Eigenspace selection criteria based upon the design goals are discussed, 
and for the study case it would appear that a desirable eigenspace can be obtained. 

In addition, the importance of state-space selection is noted along with problems 
with reduced-order measurement feedback. Since the full-state control laws may be 
implemented with dynamic compensation (state estimation), the use of reduced-order 
measurement feedback is less desirable. This is especially true since no change 
in the transient response from the pilot's input results if state estimation is 
used appropriately. The potential is also noted for high actuator bandwidth 
requirements if the linear quadratic synthesis approach is utilized. Even with 
the actuator pole location selected, a problem with unmodeled modes is noted due 
to high bandwidth. 

Some suggestions for future research include investigating how to choose an 
eigenspace that will achieve certain desired dynamics and stability robustness, 
determining how the choice of measurements effects synthesis results, and exploring 
how the phase relationships between desired eigenvector elements effects the 
synthesis results. 
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